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Summary of Progress 


In this report, we will concentrate on two aspects of our work for NASA: the construction 
of multi-dimensional phase modulation trellis codes and a performance analysis of these 
codes. Appendix A contains the final version of a paper entitled “Trellis Coded Multi- 
Dimensional Phase Modulation” [1], earlier versions of which were included in our June 1988 
and February 1987 reports. This paper, which will appear shortly in the IEEE Transactions 
on Information Theory , now contains a complete list of all the best trellis codes for use with 
phase modulation. LxMPSK signal constellations are included for M = 4, 8, and 16 and L = 
1, 2, 3, and 4. Spectral efficiencies range from 1 bit/channel symbol (equivalent to rate 1/2 
coded QPSK) to 3.75 bits/channel symbol (equivalent to rate 15/16 coded 16-PSK). The 
parity check polynomials, rotational invariance properties, free distance, path multiplicities, 
and coding gains are given for all codes. These codes are considered to be the best candidates 
for implementation of a high speed decoder for satellite transmission. In our next report, 
we will discuss the design of a hardware decoder for one of these codes, viz., the 16-state 
3x8-PSK code with free distance 4.0 and coding gain 3.57 dB. This work is being conducted 
by Mr. Steven S. Pietrobon, a Ph.D. student supported by the grant. 

Appendix B contains an exhaustive simulation study of the multi-dimensional phase 
modulation trellis codes. This study was motivated by the fact that coding gains quoted 
for almost all codes found in the literature are in fact only asymptotic coding gains, i.e. , 
the coding gain at very high SNR or very low BER. These asymptotic coding gains can be 
obtained directly from a knowledge of the free distance of the code. On the other hand, 
real coding gains at BER’s in the range of 10 -2 — 10 -6 , where these codes are most likely 
to operate in a concatenated system, must be obtained by simulation. Our study goes far 
beyond anything previously published on this subject and forms the basis for a paper recently 
presented at the Allerton Conference on Communication, Control, and Computing [2]. This 
study was conducted by Mr. Lance C. Perez, a Ph.D. student supported by the grant. 

Much other work has also been continuing since our last report. We will not discuss these 
in any detail here, but important results will be included in future reports. Other publications 
and presentations supported by the grant are listed in the references [3-9]. Copies of these as 
well as reprints of papers previously published with grant support are being sent to NASA 
under separate cover. 
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Abstract 

In this paper, trellis coded multi-dimensional MPSK modulation is investigated. A 2L- 
dimensional MPSK (LxMPSK) signal set is obtained by forming the Cartesian product of 
L 2-dimensional MPSK signal sets. A systematic approach to partitioning LxMPSK signal 
sets is used which is based on block coding. An encoder system approach is developed which 
incorporates the design of a differential precoder, a systematic convolutional encoder, and 
a signal set mapper. Trellis coded Lx4PSK, LxSPSK, and Lxl6PSK modulation schemes 
are found for 1 < L < 4 and a variety of code rates and decoder complexities, many of which 
are fully transparent to discrete phase rotations of the signal set. The new codes achieve 
asymptotic coding gains up to 5.85 dB. 

1 This work was supported in part by NASA Grant NAG5-557 and OTC (Australia) RAD Programme 
No. 4. 
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1 Introduction 


Since the publication of the paper by Ungerboeck [1], Trellis Coded Modulation (TCM) has 
become a very active research area [2-13]. The basic idea of TCM is that by trellis coding 
onto an expanded signal set (relative to that needed for uncoded transmission), both power 
and bandwidth efficient communication can be achieved. 

TCM can be classified into two basic types, the lattice type (e.g., M-PAM and M-QASK) 
and the constant amplitude type (e.g., MPSK). Constant amplitude modulation schemes 
have a lower power efficiency compared with lattice type modulation schemes but are more 
suitable for certain channels, e.g., satellite channels containing nonlinear amplifiers such as 
traveling wave tubes (TWT). Taylor and Chan [5] and Wilson et. al. [6] have studied the 
performance of trellis coded MPSK (TC-MPSK) modulation, in particular rate 2/3 TC-8PSK 
and rate 3/4 TC-16PSK, respectively, for various channel bandwidths and TWT operating 
points. Their results showed that TC-MPSK modulation schemes are quite robust under 
typical channel conditions. 

In any TCM design, partitioning of the signal set into subsets with increasing minimum 
intra-subset distances plays a central role. It defines the signal mapping used by the mod- 
ulator and provides a tight bound on the minimum free Euclidean distance (d/ ree ) between 
code sequences. For lattice-type TCM, Calderbank and Sloane [10] have made the important 
observation that partitioning the signal set into subsets corresponds to partitioning a lattice 
into a sublattice and its cosets. Forney [13] has developed a method, called the squaring 
construction, of constructing higher-dimensional lattices from partitioned lower-dimensional 
lattices. 

In this paper, we investigate a class of trellis coded multi-dimensional (multi-D) MPSK 
modulation schemes. Signals from a 2L-dimensional (2L-D) MPSK signal set (which we shall 
denote as Lx MPSK) are transmitted over a 2-D modulation channel by sending L consec- 
utive signals of an MPSK signal set. Therefore, the LxMPSK signal set is the Cartesian 
product of L 2-D MPSK signal sets. Trellis coded multi-D phase modulation (TC-LxMPSIv) 
provides us with a number of advantages that usually cannot be found with TC-MPSK: (i) 
flexibility in achieving a variety of fractional information rates, (ii) codes which are partially 
or totally transparent to discrete phase rotations of the signal set, (iii) suitability for use as 
inner codes in a concatenated coding system [14], due to their byte oriented nature, and (iv) 
higher decoder speeds resulting from the high rate codes used (rate k/(k + 1) with k up to 
15 for some codes). 

In Section 2, we introduce a block coding technique for partitioning LxMPSK signal sets. 
Section 3 describes how the encoder system, comprising a differential precoder, a systematic 
convolutional encoder, and a multi-D signal set mapper, is obtained for the best codes found 
in a systematic code search. The signal sets are designed such that the codes can become 
transparent to integer multiples of 360°/M rotations of the MPSK signal set. Also, due to 
the way in which they are mathematically constructed, a signal set mapper can be easily 
implemented by using basic logic gates and L bit binary adders. The systematic code search 
is based on maximizing dj Tee (and thus the asymptotic coding gain) as well as minimizing the 
number of nearest neighbors ( N j ree ) for various degrees of phase transparency. TC-Lx4PSK, 
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TC-Lx8PSK, and TC-Lxl6PSK codes for L = 1 to 4 are found. For TC-LxSPSK and TC- 
Lxl6PSK, asymptotic coding gains up to 5.85 dB compared to an uncoded system are 
obtained. The TC-Lx4PSK codes exhibit asymptotic coding gains up to 7.8 dB. Among 
the L — 1 codes listed are some new codes which have improvements in N j Tee and phase 
transparency compared to codes found previously [1, 4, 6, 15]. Viterbi decoding of TC- 
LxMPSK is also discussed, concentrating on maximum likelihood decoding of the parallel 
transitions within a code trellis. 

2 Multi-D Signal Set Partitioning 

In order to describe set partitioning we will start with the familiar partitioning of the 8PSK 
signal set. This is followed with an example of multi-D signal set partitioning using the 
2x8PSI< signal set. Generalizations will be gradually introduced, so that by the end of this 
section the reader should become thoroughly familiar with the concepts involved. 


2.1 Partitioning the 8PSK Signal Set 

In partitioning the SPSK signal set, or 1 xSPSK, we form a minimum squared subset distance 
(MSSD) chain of 8 2 = 0.5S6,<^ = 2,8\ = 4, and 8 2 = oo (assuming that the average signal 
energy is one). Figure 1 illustrates this partitioning, in which each subset is equally divided 
into two smaller subsets such that the MSSD in each smaller subset is maximized. Partitoning 
continues in this manner until we have eight subsets, each containing a single point, hence 
8\ = oo. 

2.2 Partitioning 2x8PSK 

A 2 xSPSK signal set ( L = 2) is illustrated in Figure 2. We use integers y\ to indicate the 
first SPSK point and y 2 for the second SPSK point, where j/ 1 , 1/2 £ {0,1,..., 7}. Natural 
mapping is used to map the integer y 3 into each complex valued SPSK signal, i.e., y 3 
exp[V^M?/_,7r/4], for j = 1,2. We can also represent j/ x and y 2 in binary form as the vector 
l Jj = [y],y},y% with y) e {0, 1}, and where y 3 = 4? /] -f 2 y) + y °, for j = 1,2. That is, the 
least significant bit (Isb) of y 3 corresponds to the right most bit and the most significant bit 
(msb) to the left most bit. We will use this convention throughout the paper. 

To represent a 2 xSPSK signal point we form the 2x3 binary matrix, 


yi 


\yl 

y\ 

yi' 

. y ^ . 


. yl 

y \ 

yl . 


Since there are a total of six bits used to describe a signal point, the unpartitioned signal 
set (indicated by f2°) has a total of 2 6 = 64 points. We also say that D° is at partition level 
p — 0. It can easily be seen that the MSSD at partition level p — 0 is Aq = d>o, = 0.586 (we 
use large A to indicate the MSSD’s for L > 1 and small 8 for L = 1). The next partition (at 
partition level p — 1) divides into two subsets of 32 points each. We call ft 1 the subset 
that contains the all zero element (i.e., y l = y 2 = 0). The other subset of 32 points is its 
coset, labeled fF (1). In forming these two subsets, we would like their MSSD, Af, to be 
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larger than Ag. If this were not possible, then we should find a partitioning that leads to a 
maximum reduction in the number of nearest neighbors within the smaller subsets (i.e., the 
average number of signal points that are distance A? away from any point). In principle, 
the partitioning could be carried out in this heuristic manner. 

A more efficient way of partitioning A 0 is to require the column vectors of y, i.e., y' = 
[y^yh] 7 , f°r o < i < 2, to be codewords in a block code. This representation using block 
codes is also known as multilevel coding (first described by Imai and Hirakawa [16] and later 
applied to QAM by Cusack [17]). To express this mathematically, we need to introduce some 
further notation. We define C m . as that block code which contains the column vectors y’, 
for 0 < i < 2. Thus, C mo contains the least significant bits of j/i and y 2 , C mi contains the 
middle bits of y\ and y 2} and so on. The actual value of m, indicates which block code is 
being used. For L = 2 there are only three block codes that are of interest to us: Co, which 
is the (2,2) block code with Hamming distance d 0 = 1 (and code words [00] , [01] , [ll] r , 
and [10] T ), Ci, which is the (2,1) block code with Hamming distance di — 2 (and code words 
[00] r and [11] T ), and C 2 , which is the (2,0) block code having only one code word, [00] r 
and Hamming distance d 2 — oo. 

Also, since C mj denotes a block code with 2 L ~ m ' code words, we can write that the 
partition level p is the sum of all the m,’s that produce the subset fl p , i.e., p = X^_ 0 m,. 
Since there are I = log 2 M bits needed for each MPSK point, p can range from 0 to IL (0 to 
6 in this case). A shorthand way of writing which column vectors y‘ belong to which block 
codes is fi(C m2 ,C m ,,C mo ). Thus, we can write fl° = f!(C 0 , C 0 , C 0 ). Since C 0 contains all 
possible length two binary vectors, then Q° is generated. 

To obtain the next partition (at level p = 1), we let Q 1 = f2(C 0 , C 0 , Cx). This partition 
satisfies our previous comments on partitioning. That is, there are only two code words 
in C x (reducing the number of points to 32), and C x contains the all zero code word. In 
partitioning, we also require the property that all the points in fl 1 belong to Q° (written as 
ft 1 C f 1°). For this example, since Cj C C 0 , this property is satisfied. This can be stated 
more generally as fl p+1 C fl p , for 0 < p < IL — 1. Thus, if we have two partition levels p 
and p\ and p’ = p + 1, then C m < C C m< for 0 < i < I — 1. 

The partition fl 1 is equivalent to forcing the lsb’s of ?/ x and y 2 to be either both zero or 
both one. By inspection of Figure 2 we can thus see that Aj = 26 q = 1.172. In fact, we 
can use a more general expression which gives a lower bound on the MSSD. From [18,19] we 
have, 

A p > min(^|_ 1 d m/ _ 1 , . . . , 6\d mi , 6%d mo ), (1) 

where d mi is the Hamming distance of the code C mi , for 0 < i < 7— 1. From (1), we obtain 
for 2x8PSK, 

Aj > rnin(4d mj ,2d m] ,0.586d mo ). (2) 

For p = 0 and 1, we can see that (2) is satisfied with equality. In fact, due to the symmetry of 
the 8PSK signal set, (2) is an equality for all values of p. It can be seen that in partitioning 
into fF and its coset H 1 (1), we could have formed Q(Co, Ci, Co) or f!(Ci, Co, C 0 ) instead 
of fl(C 0 , C 0 , Cj). However, both these other partitions have Aj = 0.586, and are therefore 
not good partitions, since we want Aj to be as large as possible. This is because dj ree can 
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be lower bounded by 2A 2 for many trellis codes [1]. 

Ignoring for the moment how the cosets are formed, we can partition fl 1 into ft 2 and 
its coset A 2 (2), and so on. (The value within the brackets of the coset will be explained in 
Section 2.3.) Every time we partition, we want to make A 2 as large as possible. To do this 
we use the following rule. The C m> that we partition (into C mi +i) from level p to level p-\- 1 
should be the i corresponding to the smallest 8 2 d m , at partition level p. If there are two or 
more 8 2 d mi that have the smallest value, we choose the one with the smallest i. 

Note that once C m . has been partitioned to C 2 (or C /_ x in general), then that particular 
block code cannot be further partitioned (since it contains only one code word). Table 1 
illustrates the partitioning of the 2x8PSK signal set. The arrows show which C m , s are 
being partitioned as p is increased. The values of A 2 are also shown. Note that at p = 3, 
we have 8 2 d mi — 4 for both i — 1 and 2. As indicated by the above rule, i — 1 is chosen 
to be partitioned to form 0 4 . Even though A 2 = A3 = 4, partition level 4 is still useful for 
coding since the number of nearest neighbors for fl 4 is less than for H 3 . This will become 
more apparent when the actual codes are found. 

The above rule usually works quite well. For L = 3, though, some of the best partitions 
do not follow this rule. Instead, we can allow a A 2 to be smaller than the rule proposes, in 
order to obtain a larger A 2 , for some p' > p , than is possible by following the rule. 


2.3 Formation of Cosets 

Now consider partition level p — 1. We have shown that there are two subsets, namely fl 1 
and its coset Q 1 (l). To obtain 0 1 (1), we must look at how coset codes are derived from 
block codes. Recall that C x is the (2,1) block code with Hamming distance d\ = 2. The 
coset of this code, Ci(l), is formed by adding modulo-2 a non-zero code word that belongs 
to C 0 , but does not belong to Ci (called the generator r°), to all the code words in Ci. 
We illustrate this with an example. C 0 has code words [0 0] r , [0 l] r , [1 0] 7 ’, and [1 1] T 
(remember that these code words correspond to column vectors of y) and Ci has code words 
[0 0] T and [1 1] T . Therefore, the generator r° could equal [0 l] 7 or [1 0] T . We arbitrarily 
choose r° = [0 1] T . Thus, C x (l) = C x © r° = {[0 l] r ,[l 0] T }. (In this paper the symbol 
© will be used to denote modulo-2 (exclusive-OR) arithmetic and + to denote integer or 
modulo-M, M > 2, arithmetic.) Note that if r° = [1 0] T , the same coset vectors would 
have been found, except that they would have been in a different order. Also note that the 
Hamming distance between codewords in C x (l) is equal to d\. 

We can also write a general expression for the cosets at partition level p = 1 as 

Ci(C°) = Ci ©C°t°, (3) 

where C° 6 {0, 1}. Thus, when £° = 0, we obtain C x (0) = C x and when C," 0 = 1, we obtain 
the coset of C X ,C X (1 ). In a similar way we can divide C x into C 2 and its coset C 2 (2), and 
C x (l) into cosets C 2 (l) and C 2 (3). Figure 3 gives an illustration of this partition. For the 
second generator, there is only one choice, i.e., r 1 = [1 1] T . The general expression for the 
cosets at partition level p — 2 becomes 

c 2 (2C' + C°) = c.ecVecV 
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= C‘[J]®C°[J], (4) 

where C2 is the all zero vector and £ {0,1}, for 0 < m < 1. We also note that 

C2 C Ci C Co and that r m £ C m , but that r m £ C m +i, for 0 < rrt < 1 . 

Since we have shown how the cosets of C m are formed, we can now show how the cosets 
of ft p are formed. We start with the simplest case, the single coset of ft 1 namely In 

the same way as the block codes are partitioned, we must find a 2 x 3 matrix that belongs 
to ft 0 but does not belong to ft 1 . This is called the generator of ft 1 and is labeled t°. Since 
C mo is partitioned in going from 0 ° to ft 1 , this implies that t° = [0,0, r°], where 0 is the all 
zero vector [0 0] r , he., 



An alternate notation for t° (using the symbol t°), is to treat t° as if it represented two 
integer values, y\ and 1/2- Thus, t° in integer form is t° = [0 1 ] T . 

To form the coset ft^l), all that is required is to add t° modulo-2 to all the signal points 
in ft 1 . We write this as 

ftV) -ft 1 ®^ 0 , (5) 

where z° £ {0,1} indicates which of the two subsets is being selected. We can see that 

in coset ft 1 (l), the Isb’s of j/i and 2/2 are either 0 and 1 or 1 and 0 , respectively. Thus this 
coset has the same MSSD as ft 1 , i.e., A 2 = 1 . 172 . Alternately, t° can be added modulo-M 
(modulo-8 in this case) to the signal points in ft 1 . With modulo-8 arithmetic, the Isb’s of y 1 
and 3/2 are still added modulo-2, but the Isb’s now produce carries which affect the middle 
and most significant bits. This is denoted as 

ft 1 ^ 0 ) = ft 1 -f z°t° (mod 8). (6) 

For example, a signal y = [1 3 ] 7 (where y = [yiJ/2] T ) in ft 1 becomes [1 2 r ] with modulo-2 
addition of t° to y or [1 4 ] r with modulo-8 addition of t° to y. Using either type of arithmetic, 
we still obtain the required partition, although the ordering of signal points within each coset 
is different. In constructing rotationally invariant trellis codes, we will find that there is a 
distinct advantage to using modulo-M arithmetic over modulo-2 arithmetic. 

Continuing with the set partitioning, it should be obvious that the next generator is 
t l — [1 1 ] T . From Table 1, we see that t l corresponds to the generator of C x . The expression 
for the cosets of ft 2 is 

ft 2 ( 2 ^ 1 + z°) = ft 2 + z 1 | + z° ^ (mod 8), ( 7 ) 

where 2* £ { 0 , 1 }, for 0 < i < 1 . For partition level p = 3 , we choose t 2 = [0 2 ] T , with 
z 2 £ { 0 , 1} used to select t 2 . Continuing in the same way, we can partition the signal set 
until we obtain only a single ( 4 -D) signal point. Thus we can form the equation (using the 
generators from Table 1 ). — 

z) = 2/1 1 = n% 2 ) 

_ 2/2 I 
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j +2° ° (mod8), (S) 

where 2 = ELo 2 ’^ » with 2 ' € {0,1}, for 0 < * < 5, and y(z) gives the integer repre- 
sentations of the two 8PSK signal points. The signal set mapping given by z can now be 
directly used by a convolutional encoder. Since y\ and ?/2 can be described in terms of 2 , the 
signal set mapper can be implemented using simple logic circuits (exclusive-OR circuits for 
modulo-2 addition and binary adders for modulo-M addition). Alternatively, since z can be 
represented with only six bits, one can use a small ROM. Figure 4 illustrates two possible 
signal set mappers for 2x8PSK. Figure 4(a) shows a mapper using modulo-2 arithmetic, and 
Figure 4(b) shows a mapper using modulo-8 arithmetic. 

In general, we can write (8) as 

2/i il - 1 

y(z) = : = (l IL (z) = £ z't\ ( 9 ) 

1=0 

L Vl J 

where 2 = Ef=< 0 1 2'V, with z* G {0,1}, for 0 < i < IL - 1. The addition in (9) is not 
specified, but may be modulo-2 (using the binary matrix generators), modulo-M (using the 
integer generators), or a combination of modulo-2 and modulo-M. Figure 5 illustrates the 
partitioning of fl° into H 3 and its cosets fl 3 (4 z 2 + 2 2 1 + z° j for the 2x8PSK signal set using 
modulo-8 addition. 

2.4 Partitioning 3xMPSK and 4xMPSK Signal Sets 

In a similar fashion to 2x8PSK, to partition LxSPSK (for L > 2) requires the partitioning 
of length L > 2 block codes. We again look for partitions that have an increasing Hamming 
distance. For L = 3, there are two partitions that are interesting. 

The first partition has Hamming distances d 0 = l,d} = 2, d\ = 2, and c/ 3 = 00 . These 
Hamming distances correspond to the (3,3), (3,2), (3,1), and (3,0) block codes Co, C}, C 2 , 
and C 3 , respectively, where C 3 C C) C Cj C C 0 . Table 2(a) gives the three generators, 
Tj 0 , r-j 1 , and T 3 , that were chosen, along with the Hamming distances (d m ) and the number of 
nearest neighbors (N m ) at each partition level m. The choice was not completely arbitrary, 
since one of the generators must be the all ones vector (which in this case is r 3 ). The reason 
for this will be explained in Section 3. 

It is interesting to note that the generator matrix for these block codes can be formed from 
the generators. In general, a generator matrix G m for an (L,L — m) block code C m , for 0 < 
m < L — 1, can be formed from the generators r m to r L-1 , i.e., G m = [r m , r m+1 , . . . ,r i,_1 ] T . 
For example, for the L — 3 block codes given in Table 2(a), 

1U 1 r 1 1 0 1 

110 ,G{= LlU , G, = [Oil]. 

0 1 1 J L J 

For the other L = 3 partition, we have d 0 = l,c/ 2 = l,d\ — 3, and d 3 = 00 . These 
distances correspond to block codes Co,C 2 ,C 2 , and C 3 , where C 3 C C 3 C Cj C Cq. Table 
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2(b) shows the generators for these codes. Note that is the all ones vector in this case. 
The advantage of this partition is that d\ = 3 is larger than d\ = 2. However, d\ = 1 is less 
than d\ =2. 

The partitions of 3x8PSK that will be useful for trellis coding are given in Table 3. 
Table 3(a) corresponds to the first partition where we try to maximize at each partition 
level. In Tables 3(b) and 3(c), the second set of block codes are used to increase Aj to 1.757 
while Aj decreases to 0.586. In Table 3(c), Ag increases to 6.0 and A4 decreases to 2.0. 
Note how A| = 6.0 is obtained in Table 3(c). At p = 4 we have A^ = min(4.0, 2.0, 00) and 
at the next partition level, A? = min(4.0, 6.0, 00) = 4.0. Now C m2 is partitioned to give 
A| = min(8.0, 6.0, 00) = 6.0. In the next level, we partition C m , to obtain A| = 8.0. In 
Section 3, the reasons why these latter two partitions are used will be seen more clearly. 

For 7 = 4 there is only one good way to partition length 4 block codes. Table 2(c) gives 
a summary of the basic parameters. Using Table 2(c), we can partition the 4x8PSK signal 
set as shown in Table 4. 

For Lx4PSIv and Lxl6PSIv we obtain from (1) that, 

Ap > rnin(4d m ,,2d mo ), (10a) 

A l > min(4d ma , 2d m2 , 0.586d m , , 0.152d mo ), (10b) 

respectively, where p = £a=d m i(I — 2 for (10a) and 7 = 4 for (10b)). In a similar fashion 
to Lx8PSK, the signal set partitions can be obtained for 7 = 2 to 4. Tables 5, 6, and 7 give 
a summary of the partitions for Lx4PSK, LxSPSK, and Lxl6PSK, respectively. 

2.5 Larger Dimensional MPSK Signal Sets and the Squaring 
Construction 

One way to obtain larger dimensional MPSK signal sets is to take an Lx MPSK signal set 
partition (with its corresponding MSSD’s relabeled as 6f, for 0 < i < IL ) and form a 277' 
dimensional MPSK signal set which we label as L' x Lx MPSK. Thus if we have a 2x8PSK 
signal set, the MSSD’s Ap,0 < p < 67', for L' x 2x 8PSK are given by 

Ap > min(8d mj , 4d m< , 4d m3 , 2d mj , 1.172d m ,,0.586d mo ), (11) 

where the d m ds are the Hamming distances of ( L L'—rrij) block codes. If L 1 = 2 we can form 
the 2x2x8FSK signal set, which is equivalent to the 4x8PSK signal set. Table 8 illustrates 
this partitioning. Note that the MSSD’s obtained are exactly the same as those found with 
the 4x8PSK partitioning given in Table 4. Figure 6 shows a block diagram of a signal set 
mapper for the partition of 2x2x8PSK. The function T\ corresponds to the mapping given 
by the generators in Table 8 and T 2 to the generators in Table 4. 

For L' = 2, the above method of obtaining larger dimensional MPSK is essentially equiv- 
alent to the squaring or two-construction described by Forney [13]. The cubing or three- 
construction corresponds to L' = 3. One can continue squaring or cubing various multi-D 
signal sets in an iterative fashion to obtain many larger dimensional signal sets. If we desire 
an Lx MPSK signal set, all that is required is to factor L to determine which constructions 
are needed. For example, if 7 = 24, we could factor this into a 2 x 2 x 2 x 3x8PSI< signal set. 
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If L is a prime number, then the appropriate length L block codes and their corresponding 
generators must be found. 

Table 9 gives the generators for L — 5 and 7. Also given are the Hamming distances and 
the number of nearest neighbors for each length L block code. Note that there are three 
different partitions for L — 5 and four different partitions for L = 7. This seems to suggest 
that the number of useful partitions increases by one for each successive prime number. 
Thus, L = 11 is expected to have five useful partitions, and so on. These partitions were 
constructed by hand and probably represent the practical limit of hand constructions. For 
L = 11 and above, an algorithmic or mathematical method is required. In forming each 
partition, we have tried to maximize the Hamming distance and minimize the number of 
nearest neighbors. For example, the type IV partition maximizes the Hamming distance 
and minimizes the number of nearest neighbors for the (7,4) block code while the type III 
partition maximizes the Hamming distance and minimizes the number of nearest neighbors 
for the (7,3) and (7,2) block codes. 

For larger dimensions, these methods may produce block codes which do not have the 
largest possible minimum distance. For example, the largest Hamming distance that can be 
obtained for the (24,12) coset code is six. However, the (24,12) Golay code has a Hamming 
distance of eight. For L = 2, 3, and 4, the block codes are relatively simple. Thus, we are 
fairly certain that the best partitions for these LxMPSK signal sets have been found. 

3 Trellis Coded Multi-D MPSK Design 

This section describes how convolutional codes are constructed for the LxMPSK signal sets 
described previously. We first show how to construct signal sets which have good phase 
rotation properties. Following this, a method used to find good convolutional codes based 
on the parity check equations is presented. 

3.1 Construction of Signal Sets 

Equation (9) can be used to describe a signal point in an LxMPSK signal set. The number 
of bits z J needed to describe each signal point is IL. If the lsb is used for coding, we can 
form a rate (IL — 1 )fIL code. A more convenient measure of rate is to use the average 
number of information bits transmitted during each 2-D signal period (T). This is called the 
effective rate of the code, R e /f = ( IL — 1 )/L (bit /T). The unit bit/s/Hz can also be used 
(for the actual bandwidth efficiency), but this assumes that perfect Nyquist filtering is used 
in the receive and transmit filters. Since this is not the case in many practical systems, we 
make a distinction between the units bit /T and bit/s/Hz. 

Other rates can be achieved by setting the q lsb’s of the mapping to 0. We do this to 
insure that the MSSD’s are as large as possible, so that the best codes can be found. In this 
case (9) can be rewritten as 

IL - 1 

= £ z'-V, (12) 

i=i 
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for 0 < z < 2 IL ~ q ~ l — 1,0 < q < L — 1 , and where y q {z) represents a point 2 in an 
LxMPSK signal set such that the first q bits of (9) are 0. As before, we do not restrict the 
type of addition that is used. We now let z = [z IL ~ q ~ l , . . . , z 1 , z°], where z is the binary 
representation of z, and the lsb of z is always the coding bit. This notation insures that the 
parity check equations of a convolutional code can always be expressed in terms of the lsb’s 
of z without depending on the type of signal set used or its partitioning. From (12), codes 
with effective rates R e ff = (I L — q — l)/L can be formed. An upper limit of q = L — 1 
is set because for q > L the signal set is partitioned such that d mo = oo, i.e., an M/2 1 
-PSK, for j > 1, signal set is being used (one exception is the 4x8PSK signal set (Table 4) 
where d mo = 4 for q = L). The MSSD’s range from A 2 to Af L and the uncoded minimum 
squared Euclidean distance (MSED) is A 2 +1 , since uncoded transmission uses only half as 
many signals as coded transmission. 


Example 3.1 

We can form a rate 4/5 code with an effective rate of 2.0 bit/T from a 2xSPSK ( L = 
2,7 = 3) signal set with q = 1. Then 


y\z) = z 4 


4 

4 


+ 2“ 


0 

4 


+ 


2 

2 


+ Z 1 


+ z“ 


(mod 8). 


The uncoded MSED is A| = 2.0, which is the same as uncoded 4PSK. 


3.2 Effect of a 360 °/M Phase Rotation on a Multi-D MPSK Sig- 
nal Set 


Using modulo-M arithmetic in (12), multi-D signal sets can be constructed such that there 
are at most I bits in z affected by a signal set rotation of ^ = 360°/A7. For 4PSK, 8PSK, 
and 16PSK, this corresponds to rotations of 90°, 45°, and 22.5°, respectively. Initially, we 
consider all possible mapped bits, i.e., <7 = 0. 

Consider that a lxMPSK signal set has been rotated by VD. Since we are using natural 
mapping, the integer representation of the rotated signal point is y r = y + 1 (mod M), where 
y is the integer representation of the signal point before rotation. If y is in binary notation, 
then * 


y? 

= y° © 1 = y°, 

(13a) 

Vr 

= y 1 0 y°, 

(13b) 

Vr 

= y 2 e.y°-y\ 

(13c) 


If there are 7 = log 2 M bits in a signal set, then we see from (13) that all I bits are 
affected by a phase rotation of 'P. 

Consider the 2x8PSK signal set, with the mapping given by (8). The phase rotation 
equations of this mapping can be determined as follows. From (8), the signal outputs can 
be written in terms of z as 


V i 
1/2 


(4z 5 + 2z 3 + z 1 ) | 


+ (4z 4 + 2z 2 + z°) 


0 

1 


(mod 8). 


(14) 


10 



be written in terms of z as 


y i 
y 2 


(4z 5 + 2z 3 + z 1 ) j 


+ (4z 4 + 2z 2 + z°) 


(mod S). 


(14) 


After a 45° phase rotation we have y^ T = yj + l(mod 8), for j = 1,2. From (14), we can 
form the following phase rotation equations, 


V\ ,T 
y 2,r 


(4z 5 + 2z 3 + z 1 + 1) ^ 


+ (4z 4 + 2z 2 + z°) 


(mod 8). 


Note that a 1 is added to the term whose coset is [1 1] T . Hence this term “absorbs” the 
effect of the phase rotation, leaving the remaining term unaffected. As can be seen, bits 
z 5 ,z 3 , and z 1 are affected in a manner similar to y 2 ,y 4 , and y° in (13), and bits z 4 ,z 2 , and 
z° are unaffected by the phase rotation. Thus, we can form the phase rotation equations 



' v r 





(15) 


If the signal set had been constructed using modulo-2 addition (instead of modulo-8), only 
z° would have remained unchanged by a 45° phase rotation. Using general notation, we can 
express (14) as 


’ 2/i " 



' 1 ' 


= (2 /_1 z p/ - 1 + • 

• + 2z Pl + z *>) 

; 

. VL . 



l 


+ 2 1 1 } + ■•• +2 {$q} + {$<>} (mod A7), (16) 

where pj, for 0 < j < 7 — 1, corresponds to those partition levels where t p equals the vector 
[2 j ,2 j , . . . ,2- 7 ] t . The term g, for 0 < j < I — 1, corresponds to those remaining terms that 
have at least one (but not all) component in t p with value 2 3 . For (14) we would have 
p 0 = l,pj = 3, and p 2 = 5. These values of pj are given for all the signal set partitions 
shown in Tables 5-7. We can now write the phase rotation equations as 


Z p ° = ©l,z* =Z P1 ©z 


P 0 


,P2 -P2 


® z 


P 0 




(17) 


and for all other partition levels, z p = z p . 

For L = 2, there is only one term in each gj. However, for L > 3, there are two or more 
terms in each g r Since the terms in g : do not contribute to the phase rotational properties 
of the signal mapping, these terms can be added modulo-2 before being added modulo-M to 
the other terms. This is best illustrated with an example. For the 3x8PSK (I) signal set in 
Table 3(a), we have the following mapping equation: 


’ J/i ' 
2/2 

= z 8 

' 0 ' 
4 

+ z T 

i 

1 

+ 2 6 


+ 2 5 

' 0 ‘ 
2 

+ 2 4 

" 2 " 
2 

. 2/3 . 


_ 4 _ 


o J 


4 J 


. 2 . 


_ 0 _ 
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+ 2 " 


1 

to to 
1 

+ z 2 

■ 0 ■ 
1 

+ z 1 

‘ 1 ‘ 

1 

+ z° 

■ i ■ 
i 

L 2 J 


. 1 . 
r 

1 1 

. 0 . 

f 

r o i 

. i . 


= (4 z 6 + 2z 3 + z°) 


+ 4 


~8 


©2 7 


1 

1 

0 


+ 2 < 



■ 0 ' 


' 1 ‘ 

1 

f 

' 0 ■ 


' 1 ' 

] 


1 

© z 4 

1 


> + lz 2 

1 

© Z 1 

1 

l 

1 

k 

. 1 _ 


_ 0 . 

J 

1 

_ 1 _ 


. 0 . 

J 


(mod 8) 


= ( 4z 6 + 2z 3 + z ° ) 





z 4 


■ z 1 ' 

+ 4 

z 8 © z T 

z 8 

+ 2 

z 5 © z 4 
z 5 

+ 

z 2 © z 1 
z 2 


(mod 8). 


The reason for this combination of modulo-2 and modulo-M arithmetic is that it reduces 
the number of logic circuits required in a signal set mapper. For small IL , it may be 
simpler to use ROM’s for signal set mapping, but for large IL this dual addition becomes 
preferable. Figure 7 gives a block diagram of the three 3x8PSK signal set mappers and 
Figure 8 illustrates the mapper for 4x8PSK. This combination of modulo-2 and modulo-M 
addition has no effect on the MSSD’s (at least for L < 4). In a similar manner, we can also 
obtain the signal set mappers for Lx4PSK and Lxl6PSK. 

Due to the phase rotational properties and simplified hardware that the combined modulo- 
2 and modulo-M mapping allows, these are the signal sets that are used to find all the trellis 
codes in this paper. 

We have shown that for q = 0, the bits that are affected by a phase rotation of 'f r are z Pj , 
for 0 < j <1—1. For q > 0 the bits that are affected are z p >~ q , for 0 < j < I — 1. However, 
depending on the signal set, pj — q for some j may be less than zero. If this is true, the 
minimum phase transparency is 2 d 'ty , where d' is the number of terms Pj — q that are less 
than zero, and the number of bits that are affected by a 2 d ^ phase rotation is s' — I — d' . 
For example, the 3x8PSI\ signal set in Table 3(a) has po = 0,pi = 3, and p 2 = 6. Thus 
if q = 1, then po — q — —1, which is less than zero, implying that d' — 1, and thus only 
s' = / — d' — 2 bits are affected by a 2'£ = 90° phase rotation. (A phase rotation of 'fr = 45° 
of this signal set produces its coset.) 

Fortunately, for the codes and signal sets considered in this paper, the above complication 
does not occur. This is partly due to the fact that for many signal sets with q = 0, the first 
L — 1 lsb’s are not affected by a phase rotation of \&. Since we consider only signal sets with 
0 < — l,<i / = 0in these cases. For those signal sets where this is not true (e.g., in some 

3xMPSK signal sets), it has been found that the convolutional codes produced are inferior 
(in either df ree or number of nearest neighbors) to an alternative signal set with d' = 0. 

When a signal set is combined with a convolutional encoder we must consider the effect 
of rotating coded sequences. A similar result to above is obtained so that, depending on the 
code and the signal set, the signal set can be rotated in multiples of 2 d $ and still produce 
valid code sequences (where d defines the degree of transparency). The actual determination 
of d is described in Section 3.4. The number of bits that are affected by a 2 d phase rotation 
is s = / — d. 
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For 0 < q < L — 1, the actual bits that are affected by a phase rotation of $ are z b} , where 
bj = pj - q, for 0 < j < I - 1. More generally, the bits that are affected by a phase rotation 
of 2 rf 'i r are z C] , where Cj = Pj+d. ~ Qi for 0 < j < s — 1. These two separate notations ( bj and 
Cj ) are used because the determination of cl depends on bj, as will be shown in Section 3.4. 

3.3 The General Encoder System 

From the above information we can now construct a suitable encoder system, as illustrated 
in Figure 9. The general encoder system consists of five sections. These sections are the 
differential precoder, the binary convolutional encoder, the multi-D signal set mapper, the 
parallel to serial converter, and the 2-D signal set mapper. The convolutional encoder is 
assumed to be in feedback systematic form, as in [1]. That is, z*(D) = x^(D) for 1 < j < k, 
where D is the delay operator and polynomial notation is used. The parity sequence, z°(D), 
will be some function of itself and the x^(D), for 1 < j < k. The parity check equation 
of an encoder describes the relationship in time of the encoded bit streams. It is a very 
useful and efficient means of describing high rate convolutional codes, since it represents the 
input/output encoder relationships in a single equation. For an R = k/(k + 1) code, the 
parity check equation is 

H\D)z h (D) © • • • © H\D)z\D ) © H°(D)z°(D) = 0(D), (18) 

where k, 1 < k < k, is the number of input sequences that are checked by the encoder, 
H j (D), for 0 < j < k, is the parity check polynomial of z-*(D), and 0(D) is the all zero 
sequence. 

Since the encoder is systematic, the differential precoder codes only those bits which are 
affected by a phase rotation. The input bits into the encoder which are precoded are denoted 
u; c ° , w Cl , . . . , w 0 ’ -1 . If c 0 = 0, we replace w° (which does not exist) by z°, as shown in Figure 
9 by the dashed line (a different precoder must then be used). For example, an encoder for a 
rate S/9 code which uses the 3x8PSIv (I) signal set given in Table 3(a) may (depending on 
the phase transparency) need this modification. This is because this signal set has bo = 0, 
and thus if the code has d = 0, then z° will need to be precoded. Figure 10 illustrates the 
two types of precoders. Note that the storage elements have a delay of LT. Figure 10(a) 
illustrates the precoder with cq > 0, where there are s inputs that are precoded. The basic 
component of the precoder is the modulo-2 s adder. For most codes this is the precoder to 
be used. For the bits that are not precoded, x' — w l , for i ^ Cj. 

Figure 10(b) shows the other case, where Co = 0 and s — 1 input bits are precoded (the 
other precoded bit is z°). The adder circuit for this case is different from Figure 10(a), i.e., 
it is not a modulo-2 s adder. The Appendix gives the equations for the differential encoder 
and decoder (for both cases) and an explanation of how these circuits work. 

At this point, we summarize the notation and indicate the limits on the parameters used 
in the search for good codes. For a rate ( IL — q — l)/(IL — q ) code, 

I = no. of bits in each 2-D signal (2 < I < 4), 

M — 2 t = no. of signal points in each 2-D signal set, 
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L = no. of 2-D signal sets (1 < L < 4), 

p = partition level of signal set (0 < p < IL), 

q = the partition level p where mapping begins (0 < q < L — 1), 

z = signal set mapping parameter (0 < z < 2 p ~ q — 1), 

k = IL — q — 1 = no. of input bits to encoder, 

k = no. of bits checked by encoder (1 < k < k ), 

vp = 360°/M = minimum phase transparency with q = 0, 

Pj = the bits z Pj affected by a ^ phase rotation with q — 0, 
d = degree of phase transparency (2'*'!', for 0 < d < I), 
s = I-d = no. of bits in z affected by a 2 d '$ phase rotation, 
cj = pj +d — q = the bits z C} affected by a 2 d ^ phase rotation. 

There are two types of systematic convolutional encoders that can be constructed. Before 
proceeding with the description of these encoders, we return to the parity check equation 
given in (18). As in [1], we define the constraint length u to be the maximum degree of all 
the parity check polynomials H^(D), for 0 < j < k. For k < j < k,H^(D) = 0, since the 
bits corresponding to these polynomials are not checked by the encoder. The parity check 
polynomials are of the form 

H j (D) = 0 © h{_ 1 D u ~ x © • • • ® h[D (B h : 0 ,l < j < 1% (19a) 

H°(D) = D 1 ' ® ® • • • ® h®D © 1. < 19b > 

If k < v, we let h 3 Q = 0, for 1 < j < k. This insures that the squared Euclidean distance 
(SED) between paths in a trellis leaving or entering a state is at least A^ +1 . Thus all codes 
in this class have a MSED between all possible non-parallel coded sequences of at least 
2Ag +1 . The parallel transitions provide an upper bound on the df ree of a code. A theoretical 
justification for constructing codes in this manner can be found in [20] where it is shown, 
using random coding arguments, that these codes have a large free MSED on the average. 

A minimal systematic encoder can be implemented from (19), since /i° = 1 [lj. The 
encoding equations are 

z j (D) = x j (D),l<j<k, (20a) 

z°{D) = H k (D)x' k (D) © • • • 0 H\D)x\D) ® {H°(D) © 1 ) 2 °(D). (20b) 

An encoder implementation using (20) is shown in Figure 11. 

For all codes with v = 1 and for some codes with v > l,k = v. For these co de s we 
cannot restrict h J 0 , for 1 < j < k. This is because k checked bits require at least k terms 
in Fr j (Z)), for 1 < j < k, that are variable. If there are not enough variables, then there 
will be some non-zero x k = [x*, . . . , x 2 , x 1 ] such that H q {D)x q = 0 (mod 2). That is, 
there will be more than 2 k ~ k parallel transitions between states in the trellis. To avoid this 
problem, when k = i/, we use (19) without any restrictons. In this case, the MSED between 
all possible non-parallel coded sequences is at least A^ + A^ +1 , since the MSED between 
paths leaving a state is A^ (since h\ e {0, 1}, for 1 < j < k) and between paths entering a 
state is A^ +1 (since hi = 0, for 1 < j < k ). 
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Hill; 


signal set mapper takes the 7 bits for each 2-D signal point and produces the required real 
and imaginary (or amplitude and phase) components for a modulator. 

Example 3.2 

In this example, we describe how to implement a particular code. The code is used with 
a 3x8PSI< signal set. Thus L = 3 and 7 = 3. We also choose q = 1, so that a 2.33 bit/T 
(rate 7/8) code is formed. The partition that is used is given in Table 3(b), from which we 
obtain p 0 = 2,p x = 3, and p 2 = 6. The code is 90° transparent, so that d = 1 and 5 = 2. 
Therefore c 0 = pi - q = 2, and c x = p 2 ~ q = 5. Thus bits w 2 and w 5 are precoded using 
a modulo-4 adder. Since Co > 0, the precoder given in Figure 10(a) is used. For this code, 
k = 2 and the parity check polynomials are H°(D ) = D 4 © D 2 © D © 1,H 1 (D) = D , and 
H 2 (D) = 7) 3 © D 2 . Excluding the parallel to serial converter and the 2-D signal mapper, 
the encoder is shown in Figure 12. This code has 16 states (i/ = 4). Note that the multi-D 
signal set mapper does not correspond exactly to Figure 7(b), since <7=1. 

3.4 Convolutional Encoder Effects on Transparency 

The convolutional encoder can affect the total transparency of the system. The method 
used to determine transparency is to examine the parity check equation and the bits that 
are affected by a phase rotation. A code is transparent if its parity check equation, after 
substituting Zj(D ) with z((D), for 0 < j < k (the rotated sequences), remains the same. 
There are normally at most 7 bits that are affected by a phase rotation, z * 0 , . . . , z bl ~ x , 
bj = pj — q, for 0 < j <7 — 1. That is, 

1 z b ° = ^©l, 

2*1 = 2 6 » ©Z* 0 , 

2* 2 = ©. 2 60 • Z h , 

Assume that the largest value of bj < k is b 0 . This implies that only one term in the parity 
check equation is affected by a phase rotation. The other bits have no effect since they are 
not checked by the encoder, i.e., bj > k for 1 < j < 7— 1. The parity check equation after a 
phase rotation of if then becomes 

H\D)z\D ) © ■ ■ • © H b °(D)[z b °(D) © 1(D)) © • • • © H°(D)z°(D) = 0(7)), 
H k (D)z\D) © • • • © H b °(D)z b °(D) © • • • © H°(D)z°(D) = E[H b °(D)}(D), (22) 

where E[H b °(D)] is the modulo-2 number of non-zero terms in 77 fc °(7)) and 1(7)) = J2%-oo 
is the all one sequence (i.e., E[H b °(D)](D) = 77 6 °(Z))1(7))). Thus if there is an even number 
of terms in H bo (D), (22) is the same as (IS). That is, the code is transparent to integer 
multiples of $ phase rotations of the signal set. However, if there is an odd number of terms 
in 77*° (7)), then E[H bo (D)\ = 1 and the coset of the convolutional code is produced. Even 
though the two equations are closely related, the codes are quite different and a decoder is 
not able to produce correctly decoded data from a ^ phase rotation of the signal set. 


(21a) 

(21b) 

(21c) 
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of terms in H b °{D), (22) is the same as (18). That is, the code is transparent to integer 
multiples of 'f phase rotations of the signal set. However, if there is an odd number of terms 
in //^(D), then E[II b °(D)] = 1 and the coset of the convolutional code is produced. Even 
though the two equations are closely related, the codes are quite different and a decoder is 
not able to produce correctly decoded data from a phase rotation of the signal set. 

Now assume that the first two terms are affected by a phase rotation, i.e., the largest 
value of bj < k is b\. The terms in the parity check polynomial H b °{D)z bt> {D)®H bi ( D)z b ' ( D ) 
now become 


[H^(D) © H b '(D)]z b °(D ) © H h {D)z bi (D ) © E[H bo (D)](D). 

In this case the parity check equation is different after a phase rotation (even if E[H b °{D)} = 
0). This means that the code is not transparent to a $ phase rotation, but it could be 
transparent to 2T or 4$ phase rotations. This is because the phase rotation equations 
reduce to 


ybo bo 7 bd - 1 _ 

T * 5 ’ ' ’ 1 * 5 

z bd = z bd © 1, z b r d+1 = z 6d+1 © z bd , . . . 

for a 2 d 'f' phase rotation, where d = 1 or 2. If there is an even number of terms in H bi (D ), 
then d— 1. This is because an even number of terms in H bl (D) cancels the effect on z bl (D) 
when the signal set is rotated by 2\E That is, the code is transparent to integer multiples of 
2$ phase rotations, but not to multiples of If there is an odd number of terms in H bi (D), 
this cancellation effect does not occur, implying that d = 2 and the phase transparency is 
4 $. 

In general, if the largest value of bj < k is bj , then d — f + E\H b > (D)]. We can then 
determine those bits z C} which are affected by a 2 d \k phase rotation, i.e., Cj = bj +( i — — g, 

for 0 < j < s — 1 , where s = I — d. 

Example 3.3 

For the code given in Example 3.2, k = 2,7 = 3, and q = 1. Thus b Q = l,&i = 2, and 
&2 = 5. Since the largest value of bj < 2 is b\, then / = 1. Therefore d = 1 + E[H bl (D)] — 
1 + E[D 3 © D 2 ] — 1. Thus the code is 90° transparent, and c 0 = 2 and C\ = 5. 

3.5 Systematic Search for Good Small Constraint Length Codes 

An approximate lower bound for the symbol error probability [1] of a multi-D code is given 

b y 


where E^/Nq is the energy per information bit to single sided noise density ratio and Q(-) is 
the complementary error function. In (23), the division by L normalizes the average number 
of errors per multi-D signal to that of a 2-D signal set. 


'.(e) 


N 


free 

~L~ 


Q 


2 N 0 


\ 


(23) 
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For each multi-D signal set considered, there are a number of code rates which can be 
achieved. As u is increased, a comprehensive code search becomes time consuming due to 
the greater complexity of each code. We have thus limited our search to v + k < 10. (The 
number of checked bits k also affects the complexity of the code search.) As indicated by 

(23) , the criteria used to find the best codes are the free MSED (d 2 ree ) and the number of 
nearest neighbors ( N f ree ). We have also included the code transparency (d) as a criteria in 
the code serach. The code search algorithm that was implemented is similar to that in [1], 
but with a number of differences which include the extra criteria mentioned above. 

The actual code search involves using a rate k/(k+ 1) code. Thus two separate notations 
are used to distinguish the rate k/(k + 1) encoder and the simplified rate k/(k + 1) encoder. 
For the rate k/(k + 1) encoder, we have x„ = [ai£, . . . , X*] (the input to the encoder) and 
z n = [ 2 *, . . . , 2 *, 2 °] (the mapped bits or encoder output) at time n. Also, e n = [e£, . • • , e\, e°] 
is the modulo-2 difference between two encoder outputs z n and z' n at time n, i.e., e n = Z n ®z(,. 
Note that there are 2* +1 combinations of z n and z' n that give the same e„. For the rate 
k/(k + 1) code, we denote reduced versions of x n , z„, and e n as x n = [x„, . . . , x„], z n = 
[z'n, . . . , z\, 2 °] , and e n = [e* , . . . , e* , e°] , respectively. 

In order to find d 2 ree for a particular code, the squared Euclidean weights (SEW) w 2 (e n ) 
are used. As defined in [1], ic 2 (e n ) is the MSED between all combinations of a(z n ) and a(z'J 
such that e n = z n ® z' n and a( z n ) is the actual LxMPSK signal point. This can be defined 
as 

ta 2 (e„) = min d 2 [a(z„), a(z n ® e n )], (24) 

all Z n 

where d 2 [a(z n ),a(z' n )] is the SED between a(z n ) and a(z' n ). One can then use the all zero 
path as a reference to find d 2 ree in a code search, i.e., 

d 2 free = min ■)> ( 25 ) 

n 

where the minimization is over all allowable code sequences with the exception of the all- 
zero sequence. We can use (25) to find d 2 Tee provided that the minimization of (24) does not 
depend on 2 °, as shown by Ungerboeck [1]. 

Although the minimization of (24) does not depend on 2 ° for lxMPSK signal sets, it can- 
not be assumed that this also applies to LxMPSK for L > 2. By expressing d 2 [a(z n ), a(z n ® 
e n )] directly in terms of z n and e n , it can be shown that 3x4PSI\ (I), 3xSPSK (I and II), 
and 3xl6PSK (I, II, and III) all depend on 2 °. This implies that (25) becomes a lower bound 
in these cases. However, due to the large number of parallel transitions for these codes, we 
can still determine d 2 ree (and Nf ree ) using a slightly modified version of (25). 

Since there are 2 fc+1 values of e„, there are a total of 2 2fc+2 computations required to find 
all the values of io 2 (e n ). For example, a rate 11/12 code with 4x8PSK modulation requires 
nearly 17 million computations. This can be reduced by letting 2 ° = 0 (or 1) and minimizing 

(24) over all z n = [ 2 *, . . . , 2 *, 0]. This reduces the number of computations to 2 2k+l . In fact, 
it is possible to even further decrease the number of computations. Using some difficult 
algebraic manipulations, it can be shown that the L output bits 2 ^ corresponding to cosets 
t p with some components equal to 2 /_1 can all be set to zero. For example, the 4x8PSK 
signal set with q = 0 can have bits 2 ^, 2 ^, 2 ^°, and 2* 1 all set to 0 when minimizing (24). 
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This is due in part to the MPSK signals being antipodal for these values. Thus the total 
number of computations can be reduced to 2 2k ~ L+1 . 

In order to reduce the time needed to find d 2 ree , we note that the trellis is equivalent to a 
rate Jc/(k + l) code with 2 k ~ k parallel transitions. Also, there are 2* +1 different sets of parallel 
transitions. If the minimum SEW is found for each of these sets of parallel transitions, the 
code search is greatly simplified, since the search for a rate k/(k + 1) code is all that is needed 
and k is usually small. Thus, the SEW’s required for a rate k/(k + 1) code search are 

w 2 (e n ) = minw 2 (e„), (26) 

where the minimization is over all [e£, . . . , e£ +1 ]. We define the free MSED of this rate 
k/(k + 1) code as 

d)ree = min J2 ^(^n), (27) 

71 

where the minimization is over all allowable code sequences (e(D)) defined by 

e(D) = ejD© e 2 D 2 ®,- • ■ ® e N D N , 


for e 1? e/v ^ 0, and N > 2. The code sequences of length N = 1 are the parallel transitions, 
where the MSED is the MSSD of the parallel transitions. A code might have d 2 ree larger than 
the MSSD of the parallel transitions, implying that d 2 rce occurs along the parallel transitions. 
With k checked bits and a rate k/(k+ 1) code, the MSSD of the parallel transitions is A^ + ^ +1 . 
Thus we can express d 2 ree as 

d )ree = ^(J^, A J+* + l)' ( 2S ) 

The best value of k can be determined from the free MSED of the best code for the 
previous value of v. The search starts with v = 1 and k = 1, and we find the code with 
the best d 2 ree and N / ree . We then increase v by one and determine k as follows. If d 2 ree for 
the previous best code was d 2 ree , then k remains the same. This is because the limit of the 
parallel transitions (A 2 + ^ +1 ) has not yet been reached and the trellis connectivity needs to 
be reduced in order to increase d 2 ree or reduce N/ ree . If d 2 Tee for the previous best code was 
A 2 + £ +j , then k is increased by one from the previous value; otherwise, d 2 ree and N/ ree would 
remain the same. If d 2 ree = A^ +1 for the previous best code, then k can remain the same 
or increase by one. Both values of k should be tried in order to find the best code. The 
best code is then found for this value of v and k, and the above process is repeated for each 
increasing value of v. 

As can be seen from (24), there may be some values of e„ and z n for which w 2 (e n ) < 
d J [a(z n ), a(z„ © e n )]. The number of nearest neighbors for e n (denoted m(e n )) is defined 
as the average number of times that iy 2 (e n ) equals d 2 [a(z n ), a(z„ © e n )]. If w 2 (e n ) equals 
d 2 [a(z n ), a(z n ffie n )] for all values of z n , then m(e n ) = 1. For example, in naturally mapped 
SPSK it is found that for e n = [0 1 1] and [1 1 1], d 2 [a(z n ), a(z„ ® e n )] = 0.586 for four 
values of z„ and 3.414 for the other four values of z„. Thus m(e n ) = 0.5 for e n = [0 1 1] 
and [1 1 1]. For all other values of e n , it can be shown that m(e n ) = 1. Zehavi and Wolf 
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[21] give a general approach to determining the full code distance spectrum, whereas we are 

only interested in the number of nearest neighbors. 

We can state this generally as follows. Let the number of bits in z„ that are varied to 

find w 2 (e n ) be b. Then 

m(e n ) = E u (^ 2 ( e n) - d 2 [a(z n ), a(z„ © e n )]) 2 _i> , (29) 

where u(-) is the unit step function and the summation is over all the bits in z n that are 
varied to find u^ 2 (e n ). Normally b — k + 1, but this can be reduced to b = k — L for the 
reasons mentioned previously. 

For the simplified rate k/(k +1) code, m(e n ) is the sum of all the m(e n ) s for which 
w 2 (e n ) = w 2 (e n ), i.e., 

m(e n ) = E u ( w2 (*n) - ™ 2 (e„)) m ( e n), ( 30 ) 

where the summation is over all [e£, . . . , e£ +1 ]. We can think of m(e n ) as the total average 
number of nearest neighbors along each set of parallel transitions. 

The number of nearest neighbors for the MSSD A 2 + ^ +1 is 

N& = E w (N+*+i ~ W *M) m ( e „)> ( 31 ) 

where the summation is over all e n = [e£, . . . , e£ +1 , 0, . . . , 0]. The number of nearest neighbors 
for paths with SED d 2 ree can be calculated using m(e n ) as follows: 

Nfr'c = E n ™(e n ), (32) 

a— l n=l 

where N a is the length of a path a that has a SED of d 2 ree and A is the number of paths 
that have a SED of d 2 Tee . If d 2 ree occurs along the parallel transitions, Nf ree = A ; a> a nd 
we define the next nearest free SED and number of nearest neighbors as d 2 next = d 2 ree and 
N n ext — N f ree , respectively. (Note that d^ ext and iV ner< may not be the true next nearest 
paths, since there may be some closer paths occuring along the parallel transitions.) When 
there are several codes that have the same free MSED and number of nearest neighbors, the 
“next nearest” values are used in code selection. When d 2 Tee occurs along paths with SED 
d 2 ree ,Nf r ee = N free- The next nearest values in this case are not given in the code tables. If 

d 2 free ~ ^ 2 +j[-+i’ then N free = AA + N free- 

Example 3.4 

In Example 3.2 we have a k = 2,q = 1, rate 7/8 (2.33 bit/T) code with a 3x8PSK (II) 
signal set. After determining the mapping of the signal set, (24) was used to find the SEW’s 
for each signal point. Equation (26) determines the u> 2 (e„)’s that were used to find the best 
rate 2/3 codes. For these codes d 2 ree = A 2 + ^ +1 = A 2 = 4.0. Using (31) we determined that 
N free is 15 (after normalizing, there are only 5 paths per 2-D symbol). In the code search 
for the best rate 2/3 codes, there were many codes which had d 2 next = d 2 ree — 4.343. Thus 
(32) was used to determine N nex t f° r each best code. Table 10 gives the values of w 2 (e n ) and 
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m(e n ) for each e„ that were used in the code search. The best code with a transparency of 
90° was found to have N next = 24. 

In order to reduce the number of codes that must be tested in our code search algorithm, 
rejection rules were used. As in Rule 1 of [l], time reversal of the parity check polynomials 
was used to reject codes. Even though u> 2 (e„) and m(e n ) are used to find the best codes, 
Rule 2 in [1] can still be exploited, provided that iu 2 (e n ) = A 2 ^ n j +g ,, where r(e n ) is the 
number of trailing zero’s in e n . When this is not true, it may still be possible to find some 
combinations of the parity check polynomials that can be rejected (this was also implemented 
in our code search). Rule 3 in [1] was also used to eliminate codes. 

In the code search, a rate k/(k + 1) code is searched for a particular v. Before finding 
d 2 ree , the code search program checks to make sure that the code only produces sequences 
with length N > 2. If for some input x n ^ 0, the inputs to the systematic encoder are all 
zero, the state of the encoder goes from one state to the next as if a zero input had occurred. 
Thus parallel transitions will occur in the rate k/(k + 1) code, which should not have parallel 
transitions. Therefore, in the code search, codes at level i (1 < i < k) were rejected if for 
some [a;*,. . . , a: 1 ] / 0, x J H J (D ) (mod 2) = 0(D). 

Two programs were used in the code search, one for codes with v > k and the other 
for codes with v = k. For specific values of I, L, and q, y 9 (z), for 0 < z < 2 IL ~ q — 1, was 
generated using the coset representatives t p , for 0 < p < IL — 1, that are given in Tables 5, 
6, and 7. The squared Euclidean weights w 2 (e n ) were then calculated using (24) for all e n . 
Since the value of k can change with each v , w 2 (e n ) and m(e n ) were computed, if necessary, 
as the program went from the smallest to the largest v. 

The code search used the various rejection rules before the time consuming tasks of 
finding d 2 ree (using the bi-directional search algorithm [22]) and iV/ ree (using a technique 
based on the Viterbi algorithm). The rejection rules were organized so that the best codes 
for each of the two possible phase transparancies were found. The code search found those 
codes which had the largest free distance (for a particular transparency). If a code was found 
to have its free MSED equal to or greater than the previous best code, N / ree was determined 
and this code was listed if either its d 2 Tee or iV/ ree had improved over the previous best code. 

The octal code generators were then listed along with their d 2 ree , N j ree , and phase trans- 
parency d. A small list of codes was produced (for each code search) from which the best 
codes could be chosen. Every time that k is increased by one in the code search (which is 
done automatically), the program determines and lists A 2 + ^ +1 and N& for use in the code 
tables. 

The asymptotic coding gain 7 of each code compared to the uncoded case, as shown in 
the code tables, is 

1=101og Ud},Jdl)(dB), (33) 

where d 2 is the smallest MSSD of an equivalent uncoded 2-D or multi-D scheme. In nearly 
all cases, d 2 = A 2 +1 . For codes with a non-integer R e /f, no equivalent lxMPSK scheme 
exists which has the same R e jf, and so the equivalent uncoded multi-D signal set is used 
instead. For the 4x8PSK signal set with q — 3, R e fj = 2 b'lt/T. Thus, a natural comparison 
would be against uncoded 4PSK, which has d 2 = 2. (In this case, A 2 +1 = 2.343, which is 
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inconsistent with other codes that also have R e ff — 2 bit /T.) The asymptotic coding gains 
compared to uncoded (M/2)-PSK are found by adding to 7 the appropriate correction factor 

7M /2 = ioiogi 0 (jTTYjf) <>“*), < 34 ) 

as shown in the code tables. The transparency (in degrees) is also given for each code. The 
parity check polynomials are expressed in octal notation in the code tables, e.g., H°(D) = 

D 6 + D 4 + D 2 + D + 1 = (001 010 111)2 = (127) s - 

In Tables 11, 15, and 19, codes for TC-lx4PSK (rate 1/2 4PSK), TC-lx8PSK (rate 
2/3 8PSK), and TC-lxl6PSK (rate 3/4 16PSK), respectively, are presented. These tables 
give the best code for each phase transparency, which (to the best of our knowledge) have 
not been previously published. The best codes, without regard for phase transparency, were 
originally published by Odenwalder [15] for 4PSK (with the codes in non-systematic form), 
by Ungerboeck [1,4] for 8PSK, and by Wilson, et. al. [6] for 16PSK. 

Tables 12, 16, and 20 list the TC-2x4PSI< codes (rates of 1.5 and 1.0 bit/T), the TC- 
2x8PSK codes (2.5 and 2.0 bit/T), and the TC-2xl6PSK codes (3.5 and 3.0 bit/T), re- 
spectively. Tables 13, 17, and 21 list the TC-3x4PSK codes (1.6/, 1.33, and 1.0 bit/T), the 
TC-3x8PSI< codes (2.67, 2.33, and 2.0 bit/T), and the TC-3xl6PSI< codes (3.67, 3.33, and 

3.0 bit/T), respectively. Tables 14, 18, and 22 list the TC-4x4PSK codes (1.75, 1.5, 1.25, and 

1.0 bit/T), the TC-4x8PSK codes (2.75, 2.5, 2.25, and 2.0 bit/T), and the TC-4xl6PSK 
codes (3.75, 3.5, 3.25, and 3.0 bit/T), respectively. 

Equivalent R = 5/6, TC-2x8PSK (2.5 bit/T) codes with up to 16 states have been found 
independently by Lafanechere and Costello [8] and by Wilson [9], although with reduced 
phase transparency. The 2 state TC-Lx8PSI< and TC-Lxl6PSI< codes were also found by 
Divsalar and Simon [23]. 

In the code tables it can be seen that for the same complexity, there are usually two 
codes (and in some cases three codes) that are given. Note that the code with the worst 
phase transparency has a better free distance or a fewer number of nearest or next nearest 
neighbors. Thus, if phase transparency is not required, one should choose the less phase 
transparent code in order to obtain the maximum performance for a given complexity. 

3.6 Decoder Implementation 

When the Viterbi algorithm is used as the decoder, a measure of decoding complexity is 
given by 2 v+ ^ / L. This is the number of distinct transitions in the trellis diagram for any 
TCM scheme normalized to a 2-D signal set. The maximum bit rate of the decoder is kf d , 
where f d is the symbol speed of the decoder. Since k is quite large for multi-D signal sets 
(at least (7 — 1)T), high bit rates can be achieved. For example, a Viterbi decoder has been 
constructed for a rate 7/9 periodically time varying trellis code (PTVTC) with v = 4, k = 2, 
and 8PSK modulation [24]. This decoder has f d = 60 MHz and a bit rate of 140 Mbit/s. 
However, with the equivalent rate 7/8 code with 3x8PSI< modulation, the bit rate will be 
L = 3 times as fast, i.e., 420 Mbit/s. The branch metric calculator, though, will be more 
complicated due to the larger number of parallel transitions between states. Alternatively, 
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one could build a decoder operating at a 20 MHz speed and achieve the same bit rate of 
140 Mbit/s. In addition to providing decreased decoder complexity, this multi-D code has 
an asymptotic coding gain which is 0.56 dB greater and is 90 transparent, compared with a 
180 transparency for the PTVTC [25]. 

Although the decoding complexity of the Viterbi algorithm is measured in terms of 
2 l/+k /L, for multi-D schemes the complexity of subset (parallel transition) decoding must 
also be taken into account due to the large number of parallel transitions. 

The Viterbi decoder must find which of the 2 k ~ k parallel transitions is closest, in a max- 
imum likelihood sense, to the received signal. A brute force method would be to determine 
the metric for each of the 2 k ~ k paths and then find the minimum. This would involve at least 
2 k ~' k — 1 comparisons. Since there are 2 A+1 sets of parallel transitions, a total of 2 fc+1 - 2 k+l 
comparisons would be required. For large k and small k, this is an unacceptably large number 
of computations. 

Fortunately, as shown in [13] for binary lattices, it is possible to greatly reduce the 
number of computations required. In fact, the decoding scheme becomes very similar to 
Viterbi decoding, except that finite length sequences are used. 

To illustrate this we will present the decoding scheme for TC-2x8PSK parallel transitions 
with k = 2 and an efficiency of 2.5 bit/T (a rate 5/6 code). There are eight sets of parallel 
transitions, with eight paths in each set. Figure 13 shows the parallel transition decoding 
trellis for z = [000] (i.e., the three lsb’s are set to zero). In Figure 1, we use the notation 
A0 to indicate the whole SPSK signal set, which divides into B0 and B 1 (4PSK signal sets 
rotated 45° from each other). B0 divides into CO and C2 (2PSK signal sets rotated 90° from 
each other), and B 1 divides into Cl and C3. This notation is also used in [l] for partitioning 
an SPSK signal set. Each segment in Figure 13 thus represents two parallel lines. The length 
of this trellis equals the dimensionality L — 2 of the signal set. 

The path CO x CO corresponds to those four paths that have z 3 — 0 and C2 x C2 
corresponds to those four paths that have z 3 = 1, giving a total of eight paths. To decode, 
hard decisions can be made for CO and C2 for each time period, from which the values of z 4 
and z 5 can be determined. For example, say that CO x CO decodes into the points 04, with a 
metric of m 0 , and C2 x C 2 decodes into the points 66, with a metric of m ly where the metrics 
are the sum of the Euclidean distances (or log-likelihood metrics for a quantized channel) 
from the first and second received points. After comparing the two metrics, if m 0 < m ly 
then z 3 = 0 and the point 04 would give z 4 — 1 and z 5 = 0 (see Table 1). If mo > m i, 
then z 3 — 1, and the point 66 would give z 4 = 0 and z 5 = 1. This is equivalent to the 
add- compare- select (ACS) operation within a Viterbi decoder. 

To decode the other sets of parallel transitions,' the cosets formed by z°, z 1 , and z 2 can 
be added to the trellis paths CO x CO and C2 x C2 to form the required trellis. This is 
illustrated in Figure 14, where the ending state in the trellis indicates which set of parallel 
transitions is being decoded. In this example, there are a total of eight hard comparisons 
and eight ACS type comparisons. These 16 comparisons compare with the 56 comparisons 
required in a brute force approach, a 3.5 times reduction. 

The above maximum likelihood method can be applied to other codes where a Viterbi like 
decoder can be used to decode the parallel transitions. With this method, the complexity of 
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decoding the parallel transitions can approach the complexity of the rate k/(k + 1) Viterbi 
decoder. A simpler approach may be with large look-up tables using ROM’s. The ROM 
itself would output the k - k bits of the chosen path, along with the branch metric for that 
path. For the TC-2xPSK example given previously, we could use one ROM for each set of 
parallel transitions. If the ROM’s had eight bit words, then three bits could be used for the 
decision, and the remaining five bits for the branch metric. A total of eight ROM’s would 
then be required, one for each set of the parallel transitions. 

When using ROM’s, it is desirable to reduce the number of bits (b) required to represent 
each received 2-D signal point, since there are a total of bL bits required to address the ROM. 
One way to reduce b is to convert the ‘checkerboard’ (rectangular) type decision boundaries 
that result from separate quantization of the inphase (I) and quadrature (Q) components to 
‘dartboard’ (radial) type decision boundaries. For example, if four bits are used in I and Q 
for an 8PSK signal with checkerboard decision boundaries, a dartboard pattern as shown in 
Figure 15 may be used instead with a total of five bits to represent each point (a reduction 
of three bits). A ROM may be used to do the conversion, or the dartboard pattern may be 
already available as polar coordinates from a digital demodulator. 

A problem with TC-LxMPSK is the need to synchronize the decoder with the L 2-D 
symbols on each trellis branch. For q = 0, most codes are fully transparent. The decoder 
performance can then be used to find the correct synchronization with the received sequence. 
For q > 0, many codes are not fully transparent, and the decoder will need to synchronize 
to one of the 2 d L possibilities (which can be quite large for some codes). However, one can 
take advantage of the fact that not all signal points are used for q > 0. For example, the 
2x8PSK signal set with q— 1 consists of the signal sets BOxBO or BlxBl. The synchronizer 
would find the smallest distance between a received pair of points and the expected signal 
set. These distances would then be accumulated over a sufficient length of time to make a 
reliable decision on the symbol timing. 

If we let each signal point be represented by its phase (since the amplitude is constant 
for 8PSK), we can write BO = {0°, 90°, 180°, 270°}, and B 1 = {45°, 135°, 225°, 315°}. Let 
<f>n and (f>n represent the phase of the first and second received symbols, respectively. The 
synchronizer distance metric is then given by 

(®W- “ a l + SW- l^« “ £l) • 

*€{ 0 , 1 } \aeBi /?€£» / 

In the synchronized noiseless case, <5 n will equal zero. In the non-synchronized noiseless case, 
there are two possible outcomes for $ n , i.e., complete matchup (<!>„ = 0°) and only one signal 
is matched ($ n = 45°). If each possiblity is equally likely, then the average value of is 
22.5° . With noise, $ n can be accumulated over a sufficient length of symbols to take advan- 
tage of this average phase distance between the non-synchronized and synchronized cases to 
reliably determine symbol synchronization. This symbol synchronization is independent of 
the Viterbi decoder, so the decoder must only determine phase synchronization. 


23 


3.7 Discussion 


In order to make a comparison of all the codes listed, a plot of nominal coding gain 7* = 
10 log 10 d 2 j ree verses complexity (/? = log 2 (2 u+k /L) = v + k — log 2 £) for each code found is 
made. These plots are given in Figure 16 for effective rates of 1.0 (with 4PSK modulation), 
2.0 (8PSK), and 3.0 bit/T (16PSK), Figure 17 for effective rates of 1.5 (4PSK), 2.5 (8PSK), 
and 3.5 bit/T (16PSK), and Figure 18 (for the remaining rates). (Note that these graphs do 
not take into account the additional complexity due to parallel transitions.) Some one state 
(‘uncoded’) codes are included as well. These one state codes correspond to block coded 
(or multilevel) schemes that have recently become an active research area [26-30]. Although 
the multi-D one state codes have negative complexity (compared to trellis codes), they can 
achieve coding gains above 0 dB. 

Note from Figure 16 for TC-Lx8PSK, R e j f = 2.0 bit/T, and v = 1, that as L increases 
the complexity decreases and 7* increases, eventually reaching 6.0 dB for L = 4. Thus, for the 
8-D signal set, the complexity factor can be reduced by a factor of four, while maintaining 
u* , compared to the TC-lx8PSK code with v = 2. Beyond /? = 4 (and 7* = 6.0 dB), 
increases in asymptotic coding gain are achieved with the new codes that have been found. 
With L — 4, a ceiling of 7* = 9.0 dB will be reached due to the nature of the set partitioning. 
It would seem that very complex codes are required (0 > 15) if this 9.0 dB limit is to be 
exceeded. 

Figure 16 also shows the Lxl6PSK codes with effective rates of 3.0 bit/T. For small 0, 
the same effect observed for TC-Lx8PSK and 2.0 bit/T occurs. That is, 0 decreases and 
7* increases as L increases. Between (3 — 3 and 0 = 9, the L = 1 and L = 2 codes are very 
close. 

Figure 18 illustrates the wide range of performance that can be achieved with the codes 
found. One can choose from a high rate code with 3.75 bit/T (but requiring a large amount 
of power) to a low rate code with 1.25 bit/T. In choosing a code, a designer may start 
with a required R e /f in order to obtain a certain bit rate through a bandwidth constrained 
channel. A trade-off can then be made between decoder complexity and the reduction in 
SNR that can be achieved with the codes found. Simulations or theoretical calculations of 
a few selected codes may also be made in order to obtain a more realistic assesment of the 
performance available. 

Note that many codes have the same asymptotic coding gain for increasing complexity. In 
reality, these codes do increase in performance with increasing complexity due to a decrease 
in number of nearest neighbors. This is especially noticeable for low SNR where the effect 
of nearest neighbors becomes more important. 

4 Conclusions 

An efficient method of partitioning multi-dimensional MPSK signal sets has been presented 
that leads to easily implemeted multi-D signal set mappers. When these signal sets are 
combined with trellis codes to form a rate k/(k + 1) code, significant asymptotic coding 
gains in comparison to an uncoded system are achieved. These codes provide a number 
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of advantages compared to trellis codes with 2-D signal sets. Most importantly, R e ff can 
vary from I - 1 to / - (1 j L) bit/T, allowing the coding system designer a greater choice of 
data rates without sacrificing data quality. As R e jj approaches /, though, increased coding 
effort (in terms of decoder complexity) or higher SNR is required to achieve the same data 
performance. 

The analytical description of multi- D signal sets in terms of block code cosets, and the use 
of systematic convolutional encoding, has resulted in an encoder design (from the differential 
encoder to the 2-D signal set mapper) that allows many good codes to be found. This 
approach has also led to the construction of signal sets that allow codes to be transparent 
to multiples of 360 °/M phase rotations. In general, increasing phase transparency usually 
results in lower code performance, due to more nearest or next nearest neighbors or smaller 
free distance. 

Another advantage is decoder complexity. As a Viterbi decoder decodes k bits in each 
recursion of the algorithm, the large values of k of codes using multi-D signal sets allows 
very high bit rates to be achieved (compared to convolutional codes that map only into a 
2-D signal set). The large number of branch metric computations can be reduced either 
through the use of a modified Viterbi algorithm or large look up tables. A method has been 
presented that uses the redundancy in some signal sets to achieve symbol synchronization 
at the decoder for codes that are not fully transparent. 

Rate k/(k + 1) TC-LxMPSK codes also have the advantage of being useful as inner 
codes in a high rate concatenated coding system with Reed-Solomon (RS) outer codes over 
GF(2 k ). If the inner decoder makes errors, one trellis branch error will exactly match one 
symbol in the outer RS code word. It is shown in [14] that the symbol oriented nature of TC- 
LxMPSK inner codes can provide an improvement of up to 1 dB in the overall performance 
of a concatenated coding system when these codes replace bit oriented TC-lxMPSK inner 
codes of the same rate. 
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Appendix 


Differential Encoding and Decoding 

Let the bit streams that are differentially encoded be w c °(D), w Cl ( D ), . . . , u> Cs_1 ( D ). We first 
assume that Co > 0 (i.e. , the convolutional encoder output z°(D) is not affected by a phase 
rotation of 2 d '5, where d = / — s). Let 

5—1 

w(D) = £2V*(D). (A.l) 

i=0 

The differential encoder (or precoder) outputs are the bit streams x Co (D), x Cl (D ), . . . , 
x c, ~' (D) which go into the convolutional encoder. Similar to (A.l), we let 

x(D) = 3 '£2 i x Ci (D). (A.2) 

i=0 

For the noiseless channel, we let the Viterbi decoder output which goes into the differential 
decoder (or postcoder) be x r (D), and the output from the postcoder be uy(D). After a 2 < *'$ r 
phase rotation, we have from Section 3.2 that 


x r (D) — x(D) + 1(D)( m °d $)■> 


(A. 3) 


where S = 2 s and 1(D) is the all ones sequence. For the postcoder, we desire that w r {D) — 
w{D ) for all multiples of 2 d T phase rotations. This is achieved by defining the postcoder 
equation as 


w r (D) = ((S - 1 )D + l)x r (D)( mod S). 


(A.4) 


Substituting (A. 3) into (A.4) we obtain 


Wr(D) = ((5 — 1)D + l)(x(D) + 1(D)) 

(mod S) 

= ((5 - 1)D + l)x(D) + ((5 - 1)D + l)l(D) 

(mod S) 

= w(D) + (5 — l)l(D) + 1(D) 

(mod S) 

= w{D) + (S)l(D) 

(mod S) 

= w(D), 



as required. Notice that since 1(D) is defined to be 1 for all time, then D t l(D) = 1(D) for all 
i. In practical situations, the sequence added to x{D ) to form x r (D ) is not constant, and will 
change with time (e.g., random phase slips within a demodulator). This will introduce short 
error bursts in w r (D ) whenever a phase slip occurs due to the combined effect of decoding 
and postcoding. The precoder equation can be derived from (A.4) as 


x(D) = Dx{D) + w(D)( mod S). 


(A.5) 
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We shall now consider the case when c 0 = 0, i.e., z°(D ) is affected by a 2 d 'I/ phase 
rotation. In this case we redefine w(D ) to be 

w(D) = Y / 2 i - 1 w Ci (D), (A.6) 

t=i 

and x(D) to be 

x(D) = S '£2 i ~ 1 x c '{D). (A.7) 

1=1 

For this case, we have 2x r {D) + z°(D) = 2 x{D) + z°(D ) + 1(D), where x r (D) and z Q r (D) 
are the inputs to the postcoder for a noiseless channel. Thus, similar to (A. 4), the postcoder 
equation is defined to be 

2w r (D) = ((S — l)D + l)(2x r (D) + z°(D)) (mod S). (A.8) 

Rearranging (A.8), we obtain the precoder equation 

2x{D) = 2Dx{D) + 2w{D) + {D + S -l)z°(D) ( mod S). (A.9) 
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Table 1: 2x8PSK Signal Set Partition 


Partition 

n p 

Minimum squared (-a 2 \ 
s ubse t d i stance ^ p' 

Generator 

Level (p) 

(t p ) T 

0 

Q(C ,C ,C ) 

v 0 0 4 O 

min( 4 , 2 , 0 . 586) = 0.586 

[0 1 ] 

1 

Q(C ,C ,C ) 
o’ 0\l' 

min(4 , 2 , 1 .172) = 1 . 172 

[1 1 ] 

2 

Q(C .C .C J 

v 0 4,0 2 

min(4 , 2 ,oo) = 2.0 

[0 2 ] 

3 

Q(C ,C ,C ) 

0 4,1 2’ 

min(4 , 4 ,°o) =4.0 

[2 2 ] 

4 

Q(C ,C ,C ) 

v 4 0 2 2 

m i n( 4 ,«,«) =4.0 

[0 4] 

5 

q <C,.C 2 ,C 2 ) 

min (8 ,«>,«>) = 8.0 

[4 4] 

6 

n(c 2 ,c 2 ,c 2 ) 

min(oo ; oo,oo) — 00 

- 


Table 2: Binary generators for L = 
(a) L = 3 (I) 


m 

d 

m 

N 

m 

(X”) T 

0 


3 

[1 1 1] 

1 


3 


2 


1 



(c) L 


and 4. 


(b) L = 3 (II) 


m 



(tp T 





1 



[0 1 1] 

2 



[i 1 n 


= 4 


m 

..... j 

Efl 

N 

m 

(x m ) T 


D 

4 

[0 0 0 1] 

1 

H 

6 

[0 0 11] 

2 

H 

2 

[0101] 

3 

1 

1 

[111 1] 
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Table 3(a): 3x8PSK Signal Set Partition (I) 


Par t ition 
Level (p) 

Q p 

Minimum squared /a 2 ) 
s ubse t d i stance ^ p 

Generator 

(t p ) ¥ 

0 

Q(C ,C ,C ) 
v 0 0 .^ 0 ' 

min(4 ,2,0 .586) = 0.586 

[1 1 1 ] 

1 

Q(C ,C ,C l ) 
o o’^r 

min(4 ,2, 1 .172) = 1 . 172 

[1 1 0 ] 

2 

Q(C ,C ,C*) 

v 0 0 ^ 2' \ 

min(4 , 2 , 1 . 172) = 1 . 172 

[0 1 1 ] 

3 

n < c „.c„.c 3 ) 

min(4 , 2 ,»=) = 2.0 

[2 2 2 ] 

4 

n( c 0 ,c;,c 3 ) 

min(4 , 4 ,oo) =4.0 

[2 2 0 ] 

5 

n(c„,c',cj 

v 0^2 3 y 

min(4 , 4 ,°°) =4.0 

[0 2 2 ] 

6 

Q(C ,C ,C ) 

^0 3 3 

min( 4 ,»,») =4.0 

[4 4 4] 

7 

£2(C\C ,C ) 
1 3 3 

min ( 8 ,«’,«’) = 8.0 

[4 4 0] 

8 

Q(C' C ,C ) 

^2 3 3 

min (8 ,«>,«>) = 8.0 

[0 4 41 

9 

Q(C ,,C , f C ) 

v 3 ’ 3 ' 

min(oo,oo,oo) = oo 

- 


Table 3(b): 3x8PSK Signal Set Partition (II) 


Partition 
Level (p) 


Minimum squared ,* 2 , 
subset distance 

Generator 

(t p ) ¥ 

0 

Q(C ,C ,C ) 

v o’ o’ 4 o' 

min(4 , 2 , 0 .586) = 0.586 

[o o n 

1 

Q(C ,C ,C 2 ) 

0 0 4 , 1 

min(4 , 2 , 0 .586) = 0.586 

[0 1 1 ] 

2 

Q(C ,C ,C 2 ) 

v 0 ’ O ^2 

min(4 , 2 , 1 .757) = 1 . 757 

[i i i] 

3 

Q(C ,C ,C ) 
v 0 ^0 2 

min(4 , 2 ,°°) = 2.0 

[2 2 2 ] 

4 

Q(C ,C\C ) 
o 4 , i y 

min(4 ,4,w) =4.0 

[2 2 0 ] 

5 

Q(C .C' C ) 

min(4 , 4 ,°°) =4.0 

[0 2 2 ] 

6 

Q(C ,C ,C ) 

v ,0 3 3 

min(4 , °° , °°) =4.0 

[4 4 4] 

7 

Q(C‘ C_ ,C ) 

4, 1 J J 

min( 8 .oo.oo) = 8.0 

[4 4 0] 

8 

Q(C* C ,C ) 

4,2 3 3 

min ( 8 = 8.0 

[0 4 4] 

9 

^(C 3 ,c 3 ,c 3 ) 

min(<», 00 , 00 ) = oo 

- 


31 








Table 3(c): 3x8PSK Signal Set Partition (III) 


1 

Par t ition 
Level (p) 

a p 

Minimum squared /-a 2 -) 
s ubse t d i stance ' p' 

Generator 

(t p r 

0 

n < c ,To-c„> 

min(4 , 2 , 0 .586) = 0.586 

[0 0 1 ] 

1 

n(c o ,c o .£) 

min(4 , 2 , 0 .586) = 0 . 586 

[0 1 1 ] 

2 

Q(C ,C ,Ch 

v 0 0 

min(4 , 2 , 1 .757) = 1 . 757 

[1 1 1 ] 

3 

£2( C 0 ,C 0 .C 3 ) 

min(4 , 2 ,oo) = 2.0 

[0 0 2 ] 

4 

n( c 0 ,c;,c 3 ) 

min(4 , 2 ,°°) = 2.0 

[0 2 2 ] 

5 

£2(C„ ,C*,C,) 

min(4 , 6 ,«) =4.0 

[4 4 4] 

6 

Q(C[,C* C ) 

min( 8 , 6 , 0 °) = 6.0 

[2 2 2 ] 

7 

n(c;,c 3 .c 3 ) 

min ( 8 , °° , °°) = 8.0 

[4 4 0] 

8 

£2(C' C ,C ) 

min ( 8 , oo , o °) = 8.0 

[0 4 4] 

9 

« ( c ,c 3 ,c 3 ) 

min(t»,«>,oo) = °° 

- 


Table 4: 4x8PSK Signal Set Partition. 


Part ition 
Level (p) 

n p 

Minimum s quared / a 2 \ 
s ubse t d i s tance ' p' 

Generator 

(t p r 

0 

Q(C„ ,c .CJ 

0 0,0 

min(4 , 2 , 0 .586) = 0.586 

[0 0 0 1] 

1 

a(C ,C ,C ) 

0 Oil 

min(4 , 2,1. 172) = 1 . 172 

[0 0 11] 

2 

fl(C n ,C n ,C,) 

min(4 ,2,1. 172) = 1 . 172 

[0101] 

3 

WC 0 ,cA) 

min(4 , 2 , 2 . 343) = 2.0 

[0 0 0 2] 

4 

n(C 0 ,t r C) 

min(4 , 4 , 2 . 343) = 2 . 343 

[1111] 

5 

a< c o ,c,A> 

min(4 , 4 , <*> ) =4.0 

[0 0 2 2] 

6 

g(c„A.c 4 ) 

mi n( 4, 4, 00 ) =4.0 

[0 2 0 2] 

7 

Q(C„,C 3 , c 4 ) 

min(4 , 8 , «> ) =4.0 

[0 0 0 4] 

8 

n(t,.c,,c 4 ) 

min( 8 , 8 , 00 ) =8.0 

[2 2 2 2] 

9 

«(C,.t 4 .c 4 ) 

min(8,°o,°°) =8.0 

[0 0 4 4] 

10 

Q(C 2 . c 4 , c 4 ) 

min( 8 , °° , «> ) =8.0 

[0 4 0 4] 

11 

n<c 3 .c 4 ,c 4 ) 

min( 1 6 , °o , °°) =16.0 

[4 4 4 4] 

12 

n<c 4 ,c 4 .c 4 ) 

min(oo,o«> ,«>) =00 

- 
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Table 5: Summary of Lx4PSK partitions. 




2 


n 

L = 3 

(II) 

L=3 (III) 

L = 

■ 

Partition 

MS SD 

gen. 

H 

m 

MSSD 

gen. 

MSSD 

gen. 



Level (p) 


0 P J 

(4 p> 

(t p y 

(A 2 ) 

P 

(t p J 

(A 2 ) 

P 

(t P I 

< a :> 

|Q 

0 

2 

01 


in 

2 


2 

001 

2 

0001 

1 

4 

li 



2 


2 

Oil 

4 

0011 

2 

4 

02 



4 

222 

4 

002 

4 

0101 

3 

8 

22 


222 

6 

111 

4 

022 

4 

0002 

4 

- 

- 

8 


8 


6 

111 

8 

mi 

5 

- 

- 

8 


8 


12 

222 

8 

0022 

6 

- 

- 

- 

- ■ 

- 

- 

- 

- 

8 

0202 

7 

- 

- 

- 

- 

- 

- 

- 

- 

1 6 

2222 

Po P, 

1 

3 

0 

3 

3 

2 

n 

5 

4 

7 


Table 6: Summary of Lx8PSK partitions. 



L = 

2 

L = 3 

(I) 

L=3 

(II) 

L=3 (III) 

L = 

4 

Parti tion 

MSSD 

gen. 

MSSD 

gen. 

MSSD 

gen. 

MSSD 

gen. 



Level (p) 

< a :> 

<t p T 


(t p J 

(A 2 ) 

P 

(t p J 

<a p> 

(t p J 

(A p } 

| m 

0 

0.586 

01 

0.586 

111 

0.586 


0 .586 



0001 

1 

1.172 

11 

1.172 

110 

0.586 

Oil 



1.172 

0011 

2 

2 

02 

1.172 

011 

1.757 

111 

1 .757 

111 

1.172 

0101 

3 

4 

22 

2 

222 

2 

222 

2 


2 

0002 

4 

4 

04 

4 


4 

220 

2 


2 . 343 

mi 

5 

8 

44 

4 


4 

022 

4 

444 

4 

0022 

6 

- 

- 

4 

444 

4 

444 

6 

222 

4 

0202 

7 

- 

- 

8 

440 

8 


8 


4 

0004 

8 

- 

- 

8 

044 

8 

044 

8 

044 

8 

2222 

9 

- 

- 

- 

- 

- 

- 

- 

- 

8 

0044 

10 

- 

- 

- 

- 

- 

- 

- 

- 

8 

0404 

11 

- 

- 

- 

- 

- 

- 

- 

- 

1 6 

4444 

Po Pi P 2 

1 ; 

1 5 

0 ; 

6 

2 : 

6 

2 ( 

5 

4 8 11 
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Table 7: Summary of Lxl6PSK partitions. 
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Table 8: 2x2x8PSK signal set partition. 


P 


a p 


Mini mum 
subset 


s qua r ed 
d i s tance 


(A 2 ) 

p 


gen.. 

(t p ) 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


a ( C o- C o- C o' C o’ C o-C„) 

n < c „. c „. c „' c „. c „A> 

n < C o' C o' C o' C o-C„A> 

Q <<V C .' C o-<:.A- C 2> 

Q < C o- C o’ C oA.C,. C 2 > 

£2(C ,C ,C ,C ,C ,C ) 

Q(C ,C ,C ,C ,C ,C ) 
v o 0,0’ 2 2 r 

£2(C ,C ,C ,C ,C ,C ) 
0,0 1 2 2 2 

ft(C ,C ,C ,C ,C ,C ) 

n < c o-c,.c 2 .c 2 ,c 2 .c 2 ) 

n<c„.6 2 ,c 2 .c 2 .c 2> c 2 ) 

Q(C ] ,C,,C 2 ,C 2 ,C C ) 

n ( t 2 ,c 2 ,c 2 ,c 2 .c 2 ,c 2 ) 


min( 8 , 4, 4, 2, 1 . 172,0. 586) 

= 0 . 586 

min( 8 , 4 , 4 , 2 , 1 . 1 72, 1 .172) 

= 1 . 172 

min( 8 , 4 , 4 , 2 , 1 . 1 72, 00) 

= 1 . 172 

min( 8 , 4, 4, 2, 2. 343, 00) 

= 2.0 

min( 8, 4, 4, 4, 2. 343, o») 

= 2 . 343 

min( 8,4 , 4 , 4 , °° , «• ) 

= 4.0 

min( 8,4 , 4 ) 

= 4.0 

min( 8,4 , 8 , 00 , 00 , 00 ) 

= 4.0 

min( 8,8 , 8 

= 8.0 

min( 8,8 j 00 , 00 , 00 , 00 ) 

= 8.0 

min( 8, , 00 , °°, °°) 

= 8.0 

min( 1 6 , <*» , 00 , 00 , 00 , 00 ) 

= 16. 0 


min(°° j 00 , 00 , 00 , 00 , 00 ) = 00 


[0 1 ] 

[1 1] 

[0 2 ] 

[0 4] 

[2 2 ] 

[4 4] 

[0 8 ] 

[0 16] 
[8 8 ] 
[16 16] 
[0 32] 
[32 32] 



Table 9: Binary generators for L = 5 and 7. 



L = ‘ 

(I) 

L 


L 

= 5 (III) 

m 

d 

m 

N 

m 

dp T 

d 

m 

N 

m 

(t™) T 

d 

m 

N 

m 

(x-) T 

0 

1 

5 

[111111 

1 

5 

[11111] 

1 

5 

[00001] 

1 

2 

10 

[00011] 

1 

2 

[00001] 

1 

1 

[00010] 

2 

2 

4 

[00101] 

2 

2 

[00110] 

2 

3 

[00101] 

3 

2 

1 

[11000] 

3 

2 

[10101] 

2 

1 

[01001] 

4 

4 

1 


4 

1 

[01111] 

5 

1 

[11111] 




L = 

7 (I) 


L = 

7 (ii) 


L = 

= 7 (III) 


L = 

7 (IV) 

m 

a 

N 

m 

(x^) 1 

d 

m 

N 

m 



N 

m 



N 

m 

(^) T 

0 

n 

7 

[1111111] 

1 

7 



7 


1 

7 


1 

2 

2 1 

[0000011] 

1 

2 


1 

1 


1 

3 


2 

2 

9 

[0001001] 

2 

5 


2 

6 

[1111111] 

1 

1 


3 

2 

3 

[0010010] 

2 

1 


2 

2 


3 

7 


4 

2 

1 

[0001100] 

3 

3 


3 

2 


3 

3 


5 

4 

2 

[1111000] 

4 

2 


4 

1 


3 

1 


6 

6 

1 

[0111111] 

6 

1 


5 

1 


7 

1 

[1111 1 11] 
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Table 10: Squared Euclidean Weights used in the code search tor 
rate 7/8 (2.33 bit/T) codes with 3x8PSK (II) and k = 2. 


e 

n 

w 2 (e ) 

n 

m <*.> 

000 

0.0 

1 

001 

1.172 

2 

010 

1.757 

4 

Oil 

0.586 

1 

100 

2.0 

6 

101 

1.172 

2 

110 

1.757 

4 

111 

0.586 

1 


Table 11: Trellis Coded lx4PSK. 


R =1.0 bit JT, d 2 = 4.0, N = 1 (lx2PSK). 

eff U U 


V 

k 

h 1 

h° 

Inv. 

d 2 

free 

N 

free 

d 2 

next 

N 

next 

y (dB) 

1 

1 

1 

3 

360“ 

6 

l 

- 

- 

1.76 

2 

1 

2 

5 

360° 

10 

l 

- 

- 

3.98 

3 

1 

06 

13 

180“ 

12 

2 

- 

- 

4.77 


1 

04 

13 

360° 

12 

1 

- 

-- 

4.77 

4 

1 

06 

21 

180° 

12 

1 

- 

- 

4.77 


1 

10 

23 

360° 

14 

2 

- 

- 

5.44 

5 

1 

36 

45 

180" 

16 

2 

- 

- 

6.02 


1 

26 

53 

360° 

16 

1 

- 

- 

6.02 

6 

1 

042 

117 

180° 

20 

11 

- 

- 

6.99 

7 

1 

126 

235 

180° 

20 

2 

- 

- 

6.99 


1 

144 

223 

360° 

20 

1 

- 

- 

6.99 

8 

1 

262 

435 

180° 

24 

11 

- 

- 

7.78 


1 

362 

515 

360° 

24 

9 

- 

- 

7.78 

9 

1 

0644 

1123 

180° 

24 

2 

- 

- 

7.78 


1 

0712 

1047 

360° 

24 

1 

- 

- 

7.78 


y = 0 dB 
‘2 
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I 

Table 12(a): Trellis Coded 2x4PSK. 


R =1.5 bit/T, q=0, d 2 = 4, N =6 (2x4PSK). 

eff ^ U u 


V 

k 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N 

free 

d 2 

next 

N 

next 

Y (dB) 

1 

1 

- 

- 

1 

3 

180° 

4 

2 

6 

8 

0.00 

2 

2 

- 

1 

3 

5 

O 

90 

6 

6 

- 

- 

1.76 

3 

2 

- 

04 

06 

11 

O 

90 

8 

5 

- 

- 

3.01 

4 

2 

- 

10 

06 

23 

0 

90 

8 

1 

10 

16 

3.01 

5 

3 

14 

30 

02 

41 

180° 

10 

8 

- 

- 

3.98 


3 

16 

24 

06 

53 

360° 

10 

7 

- 

- 

3.98 

6 

3 

030 042 014 

103 

180° 

12 

40.25 

- 

- 

4.77 


3 

076 024 010 

157 

360° 

12 

30.75 

- 

- 

4.77 

7 

3 

044 022 

114 

211 

180° 

12 

8 

- 

- 

4.77 


Y = 1.76 dB 
‘2 


Table 12(b): Trellis Coded 2x4PSK. 


R =1.0 bit/T, q=l, d 2 = 4.0, N = 1 (lx2PSK). 

cff u u 


V 

k 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N 

free 

d 2 

next 

N 

n ext 

Y (dB) 

1 

1 

- 

1 

3 

U 

90 

8 

5 

- 

- 

3.01 

2 

1 

- 

2 

5 

O 

90 

8 

1 

12 

8 

3.01 

3 

2 

04 

02 

11 

360° 

12 

5 

- 

- 

4.77 

4 

2 

14 

06 

23 

180" 

12 

1 

- 

- 

4.77 

5 

2 

30 

16 

41 

180° 

16 

8 

- 

- 

6.02 

6 

2 

036 052 

115 

180° 

16 

1 

- 

- 

6.02 

7 

2 

044 

136 

203 

180° 

20 

6 

- 

- 

6.99 

8 

2 

1 10 226 433 

180° 

24 

33 

- 

- 

7.78 


y = 0 dB 
‘2 
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Table 13(a): Trellis Coded 3x4PSK. 


R =1.67 bit/T, q=0, d 2 = 4.0, N = 15 (3x4PSK I). 

eff n u u 


V 

k 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

f ree 

N f 

free 

d 2 

next 

N 

next 

y (dB) 

sig. 

set 

1 

1 

- 

- 

1 

3 

0 

90 

4 

7 

6 

32 

0.00 

I 

2 

2 

- 

2 

1 

5 

0 

90 

4 

3 

6 

24 

0.00 

I 


2 

- 

2 

1 

5 

360° 

4 

2 

- 

- 

0.00 

II 

3 

2 

- 

04 

02 

11 

0 

90 

4 

1 

6 

6 

0.00 

III 


2 

- 

04 

02 

11 

360° 

6 

11 

- 

- 

1.76 

II 


3 

05 

04 

02 

11 

0 

90 

4 

0.25 

- 

- 

0.00 

III 

4 

2 

- 

14 

02 

21 

180° 

6 

6 

- 

- 

1.76 

II 

3 

3 

01 

02 

06 

11 

360° 

6 

4 

- 

- 

1.76 

II 

4 

3 

10 

04 

02 

21 

0 

90 

6 

5.5 

- 

- 

1.76 

III 


3 

12 

04 

02 

21 

180° 

8 

19 

- 

- 

3.01 

I 

5 

3 

24 

14 

02 

41 

180° 

8 

7 

- 

- 

3.01 

I 

6 

3 

024 042 010 

105 

180° 

8 

3 

10 

16 

3.01 

I 


Y = 2.22 dB 
*2 
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Table 13(b): Trellis Coded 3x4PSK. 


R = 1.33 bit/T, q=l, d z = 4.0, N = 3 (3x4PSK II). 

eff 1 u u 


V 

K 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N. 

free 

d 2 

next 

N 

next 

Y (dB) 

sig. 

set 

1 

l 

- 

- 

1 

3 

90 

4 

l 

8 

4 

0.00 

III 


1 

- 

- 

1 

3 

360" 

6 

7 

- 

- 

1.76 

II 

2 

l 

- 

- 

2 

5 

360° 

6 

4 

10 

9 

1.76 

II 


2 

- 

2 

1 

5 

0 

90 

6 

2 

8 

4 

1.76 

III 


2 

- 

3 

1 

5 

180° 

8 

21 

- 

- 

3.01 

I 


2 

- 

2 

1 

5 

360° 

8 

16 

- 

- 

3.01 

II 

3 

2 

- 

04 

02 

11 

0 

90 

6 

2 

8 

1 

1.76 

III 


2 

- 

02 

06 

11 

180° 

8 

3 

12 

100 

3.01 

II 


3 

06 

04 

03 

11 

O 

90 

8 

1 

- 

- 

3.01 

III 

4 

3 

14 

04 

12 

23 

O 

90 

10 

5 

- 

- 

3.98 

III 

5 

3 

30 

04 

22 

43 

0 

90 

12 

13 

- 

- 

4.77 

III 

6 

3 

036 060 026 

103 

0 

90 

12 

2 

- 

- 

4.77 

III 

7 

3 

140 

160 062 213 

O 

90 

12 

1 

14 

5 

4.77 

III 


3 

004 

154 056 207 

180° 

12 

1 

16 

128 

4.77 

III 


Y = 1.25 dB 
*2 


Table 13(c): Trellis Coded 3x4PSK. 


R =1.00 bit/T, q=2, d 2 = 4.0, N = 1 (lx2PSK). 

eff 1 u u 


V 

k 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 


d 2 

next 

N 

next 

Y (dB) 

sig. 

set 

0 

0 

- 

- 

- 

- 

90° 

6 

4 

- 

- 

1.76 

II 

i 

1 

- 

- 

1 

3 

0 

90 

6 

2 

8 

l 

1.76 

III 


1 

- 

- 

1 

3 

180^ 

8 

3 

12 

16 

3.01 

II 

2 

2 

- 

3 

2 

5 

90 

10 

4 

- 

- 

3.98 

III 

3 

2 

- 

06 

02 

11 

90 

10 

2 

- 

- 

3.98 

III 


2 

- 

02 

06 

13 

180° 

12 

5 

- 

- 

4.77 

III 

4 

2 

- 

12 

16 

21 

O 

90 

12 

1 

14 

2 

4.77 

III 


2 

- 

04 

12 

27 

180° 

12 

1 

16 

22 

4.77 

III 


3 

10 

04 

02 

21 

180° 

14 

3 

- 

- 

5.44 

II 

5 

3 

22 

16 

04 

53 

180° 

16 

2 

- 

- 

6.02 

II 


3 

24 

14 

02 

i 

43 : 

0 i 

360 

16 

1 

- 

- 

6.02 

II 

6 

3 

070 004 022 

101 i 

180° 

18 

3 

- 

- 

6.53 

II 

7 

3 

156 024 046 213 

O 1 

180 

20 

3 

- 

- 

6.99 

II 


3 

044 014 

102 

217 

360° 

20 

2 

- 

* - 

6.99 

II 


y = 0.0 dB 
'2 
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Table 14(a): Trellis Coded 4x4PSK. 


R =1.75 bit/T, q=0, d 2 = 4.0, N = 28 (4x4PSK). 

eff u u 


V 

k 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

f ree 

N f 

free 



Y (dB) 

1 

1 

- 

- 

1 

3 

u 

90 

4 

12 

6 

64 

0.00 

2 

2 

- 

2 

1 

5 

O 

90 

4 

4 

6 

48 

0.00 

3 

3 

04 

02 

01 

11 

90 

6 

28 

- 

- 

1.76 

4 

3 

10 

04 

02 

21 

0 

90 

8 

78 

- 

- 

3.01 

5 

3 

24 

14 

02 

41 

O 

90 

8 


- 

- 

3.01 

6 

3 

050 032 004 

103 

O 

90 

8 

14 

10 

160 

3.01 


y = 2.43 dB 
‘2 


Table 14(b): Trellis Coded 4x4PSK. 





= 

1.50 bit/T, 

q=i> 

d 2 = 

U 

4.0, 

N = 6 (2x4PSK). 

U 


V 

E 

h 4 

h 3 

h 2 

h 1 

h° 

Inv. 

d f 
f ree 

N f 

free 

d 2 

next 

N 

next 

Y (dB) 

1 

1 

- 

- 

- 

1 

3 

— cr 

90 

4 

4 

8 

64 

0.00 

2 

2 

- 

- 

2 

1 

5 

0 

90 

8 

78 

- 

- 

3.01 

3 

2 

- 

- 

04 

02 

11 

0 

90 

8 

30 

- 

- 

3.01 

4 

2 

- 

- 

12 

04 

23 

0 

90 

8 

16 

12 

320 

3.01 

5 

3 

- 

14 

34 

06 

41 

o 

90 

8 

6 

12 

176 

3.01 


3 

- 

04 

14 

22 

43 

180° 

8 

6 

12 

160 

3.01 

6 

4 

014 006 056 022 

103 

O 

90 

8 

2 

12 

62 

3.01 


y = 1.76 dB 
'2 


8 
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Table 14(c): Trellis Coded 4x4PSK. 


R = 1.25 bit/T, q=2, d 2 = 4.0, N = 4 (4x4PSK). 

eff u u 


V 

k 

h 4 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N r 

free 

d 2 

next 

N 

next 

Y (dB) 

1 

1 

- 

- 

- 

1 

3 

U~ 

90 

8 

30 

- 

- 

3.01 

2 

1 

- 

- 

- 

2 

5 

0 

90 

8 

14 

12 

64 

3.01 

3 

2 

- 

- 

06 

02 

11 

o 

90 

8 

6 

12 

64 

3.01 


2 

- 

- 

02 

06 

11 

180° 

8 

6 

12 

32 

3.01 


3 

- 

01 

03 

06 

11 

o 

90 

8 

2 

12 

56 

3.01 

4 

3 

- 

10 

14 

06 

21 

0 

90 

8 

2 

12 

8 

3.01 

5 

4 

10 

04 

06 

22 

41 

0 

90 

12 

8 

- 

- 

4.77 

6 

; 

4 

024 014 006 042 

103 

0 

90 

16 

109 

- 

- 

6.02 


y = 0.97 dB 
‘2 


Table 14(d): Trellis Coded 4x4PSK. 


R =1.00 bit/T, q=3, d 2 = 4.0, N = 1 (lx2PSK). 

eff u U 


V 

k 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N f 

free 

d 2 

ne xt 

N 

next 

Y (dB) 

0 

0 

- 

- 

- 

- 

— tr 

9° 

8 

14 

- 

- 

3.01 

i 

1 

- 

- 

1 

3 

180° 

8 

6 

16 

64 

3.01 

2 

2 

- 

2 

3 

5 

0 

90 

8 

2 

16 

64 

3.01 

3 

3 

02 

04 

03 

11 

0 

90 

16 

45 

- 

- 

6.02 

4 

3 

02 

10 

06 

21 

0 

90 

16 

17 

- 

- 

6.02 

5 

3 

22 

10 

06 

41 

0 

90 

16 

5 

- 

- 

6.02 

6 

3 

010 060 036 

105 

0 

90 

16 

1 

20 

4 

6.02 


y = 0 dB 
'2 


42 



Table 15: Trellis Coded lx8PSK. 


R = 2.0 bit/T, d 2 = 2.0, N = 2 (lx4PSK). 

eff u u 


V 

E 

h 2 

h 1 

h° 

j 

Inv. 

d 2 

free 

N. 

free 

d 2 

next 

N 

n ext 

Y (dB) 

1 

l 

- 

1 

3 

180° 

2.586 

2 

- 

- 

1.12 

2 

1 

- 

2 

5 

180" 

4.0 

1 

4.586 

4 

3.01 

3 

2 

04 

02 

11 

360 U 

4.586 

2 

- 

- 

3.60 

4 

2 

14 

06 

23 

180° 

5.172 

4 

- 

- 

4.13 


2 

16 

04 

23 

360 U 

5.172 

2.25 

- 

- 

4.13 

5 

2 

14 

26 

53 

180° 

5.172 

0.25 

- 

- 

4.13 


2 

20 

10 

45 

360° 

5.757 

2 

- 

- 

4.59 

6 

2 

074 

012 

147 

180° 

6.343 

3.25 

- 

- 

5.01 

7 

2 

146 

052 

225 

180° 

6.343 

0.125 

- 


5.01 


2 

122 

054 

277 

360° 

6.586 

0.5 

- 

- 

5.18 

8 

2 

146 

210 

573 

180° 

7.515 

3 .375 

- 

- 

5.75 


2 

130 

072 

435 

360° 

7.515 

1.5 

- 

- 

5.75 


y 4 = 0 dB 
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Table 16(a): Trellis Coded 2x8PSK. 


R = 2.5 bit/T, q=0, d 2 = 1.172, N = 4 (2x8PSK). 

eff u u 


V 

k 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N r 

free 

d 2 

next 

N 

next 

y (dB) 

1 

1 

- 

- 

1 

3 

— cr 

90 

1.757 

8 

2.0 

4 

1.76 

2 

1 

- 

- 

2 

5 

0 

90 

2.0 

4 

2.929 

32 

2.32 

3 

2 

- 

04 

06 

11 

0 

45 

2.929 

16 

- 

- 

3.98 

4 

2 

- 

16 

12 

23 

45 

3.515 

56 

- 

- 

4.77 

5 

2 

- 

10 

06 

41 

0 

45 

3.515 

16 

- 

- 

4.77 

6 

2 

- 

004 030 

113 

0 

45 

4.0 

6 

4.101 

80 

5.33 


2 

- 

044 016 

107 

0 

90 

4.0 

6 

4.101 

48 

5.33 

7 

3 

110 

044 016 

317 

O 

90 

4.0 

2 

4.101 

25 

5.33 


y 4 = -1.35 dB 


Table 16(b): Trellis Coded 2x8PSK. 


R = 2.0 bit/T, q=l, d 2 = 2.0, N = 2 (lx4PSK). 

eff u u 


V 

k 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N 

free 

d 2 

next 

N 

next 

y (dB) 

1 

1 

- 

- 

1 

3 

- - cr 

45 

3.172 

8 

4.0 

6 

2.00 

2 

1 

- 

- 

2 

5 

0 

45 

4.0 

6 

5.172 

32 

3.01 

3 

2 

- 

04 

02 

11 

180" 

4.0 

2 

5.172 

16 

3.01 

4 

3 

04 

14 

02 

21 

0 

90 

5.172 

8 

- 

- 

4.13 

5 

3 

24 

14 

06 

43 

0 

90 

6.0 

6 

- 

- 

4.77 

6 

3 

012 050 004 

125 

0 

90 

6.343 

5.5 

- 

- 

5.01 

7 

3 

110 044 016 

317 

0 

90 

7.515 

25 

- 

- 

5.75 


y 4 = 0 dB 
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Table 17(a): Trellis Coded 3x8PSK. 


R = 2.67 bit/T, q=0, d 2 = 1.172, N = 12 (3x8PSK I). 

eff u u 


V 

k 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N r 

free 

d 2 

next 

N 

next 

Y (dB) 

sig. 

set 

1 

1 

- 

- 

1 

3 

0 

45 

1.172 

4 

- 

- 

0.00 

II 

2 

1 

- 

- 

2 

5 

O 

45 

1.757 

16 

- 

- 

1.76 

II 

3 

2 

- 

04 

02 

11 

0 

45 

2.0 

6 

2.343 

16 

2.32 

I 

4 

3 

14 

04 

02 

21 

0 

90 

2.343 

12 

- 

- 

3.01 

I 


3 

10 

04 

02 

21 

180° 

2.343 

8 

- 

- 

3.01 

I 

5 

3 

30 

14 

02 

53 

0 

90 

2.929 

48 

- 

- 

3.98 

I 

6 

3 

050 022 006 

103 

0 

90 

3.172 

12 

- 

- 

4.33 

I 

7 

3 

056 

112 004 225 

0 

90 

3.515 

84 

- 

- 

4.77 

I 


3 

100 050 022 255 

180° 

3.515 

76 

- 

- 

4.77 

I 


y 4 = -1.07 dB 


Table 17(b): Trellis Coded 3x8PSK. 


R = 2.33 bit/T, q=l, d 2 = 1.757, N = 8 (3x8PSK II). 

eff u u 


V 

k 

h 4 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 


d 2 

next 



sig. 

set 

1 

1 

- 

- 

- 

1 

3 

O" 

90 

2.0 


2.343 

16 

0.56 

II 

2 

2 

- 

- 

3 

1 

7 

O 

90 

2.586 

6 

- 

- 

1.68 

II 

3 

2 

- 

- 

06 

02 

11 

90 

3.515 

16 

- 

- 

3.01 



2 

- 

- 

04 

02 

11 

180° 

3.757 

24 

- 

- 

3.30 

II 

4 

3 

- 

10 

04 

06 

21 

O 

45 

3.757 

12 

- 

- 




2 

- 

- 

14 

02 

27 

O 

90 

4.0 

15 

4,343 

24 

3.57 

II 

5 

3 

- 

22 

16 

06 

41 

O 

45 

4.0 

7 

- 

- 

3.57 

III 

6 

3 

- 

010 046 060 

105 

O 

45 

4.0 

3 

4.686 

8 

3.57 

III 


4 

060 024 014 002 

101 

180° 

4.0 

2 

- 

- 

3.57 

III 


y 4 = 0.11 dB 




45 










Table 17(c): Trellis Coded 3x8PSK. 


R = 2.00 bit/T, q=2, d 2 = 2.0, N = 2 (lx4PSK). 

eff ^ u u 


V 

k 

h 4 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N f 

free 

d 2 

next 

N 

next 

Y (dB) 

sig. 

set 

1 

1 

- 

- 

- 

1 

3 

180° 

3.757 

24 

- 

- 

2.74 

II 

2 

1 

- 

- 

- 

2 

5 

180° 

4.0 

15 

5.757 

144 

3.01 

II 

3 

2 

- 

- 

04 

02 

11 

O 

45 

4.0 

7 

- 

- 

3.01 

III 

4 

2 

- 

- 

12 

04 

27 

0 

45 

4.0 

3 

5.757 

32 

3.01 

III 

5 

3 

- 

14 

24 

02 

41 

180° 

5.757 

17.5 

- 

- 

4.59 

III 


3 

- 

16 

22 

06 

53 

360° 

5.757 

17 

- 

- 

4.59 

III 

6 

3 

- 

030 042 014 

103 

180° 

6.0 

11 

- 

- 

4.77 

III 


4 

014 044 024 006 

103 

180° 

6.0 

4 

- 

- 

4.77 

II 


y 4 = 0 dB 
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Table 18(a): Trellis Coded 4x8PSK. 


R = 2.75 bit/T, q=0, d 2 = 1.172, N = 24 (4x8PSK). 

eff o u 


V 

k 

h 4 h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N, 

free 

d 2 

next 

N 

next 

Y (dB) 

1 

1 

- 

- 

1 

3 

or 

45 

1.172 

8 

1.757 

64 

0.00 

2 

2 

- 

2 

1 

5 

0 

45 

1.757 

48 

- 

- 

1.76 

3 

2 

- 

04 

02 

11 

1 

45 

2.0 

8 

2.343 

64 

2.32 

4 

3 

- 10 

04 

02 

21 

O 

45 

2.343 

40 

- 

- 

3.01 

5 

3 

- 30 

14 

02 

41 

O 

45 

2.343 

8 

2.929 

288 

3.01 

6 

4 

030 020 052 014 

101 

O 

45 

2.929 

136 

- 

- 

3.98 


y = -0.94 dB 

4 


Table 18(b): Trellis Coded 4x8PSK. 


R * 2.50 bit/T, q=l, d 2 = 1.172, N = 4 (2x8PSK). 

eff ^ u u 


V 

k 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N f 

free 

d 2 

next 

N 

next 

Y (dB) 

1 

1 

- 

- 

1 

3 

H 

2.0 

8 

2.343 


2.32 

2 

2 

- 

2 

1 

5 

IS 

2.343 

40 

- 


3.01 

3 

2 

- 

04 

02 

11 

45 

2.343 

8 

3.172 



4 

3 

14 

04 

02 

21 

45 

3.172 

16 

- 


4.33 

5 

3 

24 

14 

02 

41 

45 

3.515 

64 

- 

- 

4.77 

6 

3 

014 024 042 

103 

45 

4.0 

28 

4.686 

1088 

5.33 


y 4 = -1.35 dB 
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Table 18(c): Trellis Coded 4x8PSK. 


R = 2.25 bit/T, q=2, d 2 = 2.0, N = 8 (4x8PSK). 

eff u u 


V 

k 

h 4 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N 

free 

d 2 

next 

N 

next 

Y (dB) 

1 

1 

- 

- 

- 

1 

3 

45° 

2.343 

8 

3.172 

32 

0.69 

2 

2 

- 

- 

3 

1 

5 

0 

45 

3.172 

16 

- 

- 

2.00 

3 

2 

- 

- 

06 

02 

11 

0 

45 

4.0 

28 

4.343 

64 

3.01 


2 

- 

- 

02 

06 

11 

0 

90 o 

4.0 

28 

4.686 

64 

3.01 

4 

3 

- 

04 

06 

12 

21 

45 

4.0 

12 

4.686 

32 

3.01 

5 

4 

10 

04 

06 

22 

41 

0 

45 

4.0 

4 

4.686 

16 

3.01 


y = 0.51 dB 

4 


Table 18(d): Trellis Coded 4x8PSK. 


R = 2.00 bit/T, q=3, d 2 = 2.0, N = 2 (lx4PSK). 

eff u U 


V 

k 

□ 

s 

0 

h 1 

h° 





N 

next 

Y (dB) 

1 

1 

- 

- 

- 

1 

3 

— tr 

90 

4.0 

28 

4.686 

64 

3.01 

2 

2 

- 

- 

2 

3 

5 

0 

45 

4.0 

12 

4.686 

32 

3.01 

3 

3 

- 

02 

04 

03 

11 

45 

4.0 

4 

4.686 

16 

3.01 

4 

4 

10 

04 

02 

03 

21 

45 

4.686 

8 

- 

- 

3.70 

5 

4 

02 

10 

04 

22 

41 

45 

6.343 

16 

- 

- 

5.01 

6 

4 

034 044 016 036 

107 

0 

45 

6.686 

6 

- 

- 

5.24 


4 

044 024 014 016 

103 

0 

90 

7.029 

24 

- 

- 

5.46 


y 4 = 0 dB 
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Table 19: Trellis Coded lxl6PSK 


R = 3.0 bit/T, d 2 = 0.586, N = 2 (lx8PSK). 

eff u u 


V 

k 

h 2 

h 1 

h° 

Inv. 



EM 


Y (dB) 

1 

1 

- 

i 

3 

cr 

90 

0.738 




1.00 

2 

1 

- 

2 

5 

90 

1.324 

4 

- 

- 

3.54 

3 

1 

- 

06 

13 

O 

45 

1.476 

8 

- 

- 

4.01 


1 

- 

04 

13 

90 

1.476 

4 

- 

- 

4.01 

4 

1 

- 

06 

21 

45 

1.476 

4 

- 

- 

4.01 


1 

- 

10 

23 

90 

1.628 

4 

- 

- 

4.44 

5 

1 

- 

24 

43 

O 

45 

1.781 

8 

- 

- 

4.83 


1 

- 

10 

45 

0 

90 


8 

- 

- 

5.13 

6 

1 

- 

056 

135 

O 

45 

2.0 

2 

2.085 

16 

5.33 


1 

- 

032 

107 

O 

90 

2.0 

2 

2.085 

8 

5.33 

7 

1 

- 

126 

235 

O 

45 

2.0 

2 

2.366 

16 

5.33 

8 

2 

344 

162 

717 

O 

90 

2.085 

2.938 

- 

- 

5.51 


2 

224 

112 

527 

180° 

2.085 

1 .219 

- 


5.51 


Y g = 0 dB 
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Table 20(a): Trellis Coded 2xl6PSK. 


R = 3.5 bit fT, q=0, d 2 = 0.304, N = 4 (2xl6PSK). 

eff u u 


V 

k 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N. 

free 

d 2 

next 

N 

next 

Y (dB) 

1 

1 

- 

1 

3 

45" 

0.457 

8 

- 

- 

1.76 

2 

1 

- 

2 

5 

45 

0.586 

4 

0.761 

32 

2.84 

3 

2 

04 

06 

11 

22.5° 

0.761 

16 

- 

- 

3.98 

4 

2 

16 

12 

23 

22.5° 

0.913 

56 

- 

- 

4.77 

5 

2 

10 

06 

41 

22.5° 

0.913 

16 

- 

- 

4.77 

6 

2 

004 030 

113 

22.5° 

1.066 

80 

- 

- 

5.44 


2 

044 016 

107 

0 

45 

1.066 

48 

- 

- 

5.44 

7 

2 

074 

132 

217 

22.5° 

1.172 

4 

1.218 

228 

5.85 


y = -2.17 dB 
'8 


Table 20(b): Trellis Coded 2xl6PSK. 


R = 3.0 bit/T, q=l, d 2 = 0.586, N = 2 (lx8PSK) 

eff 1 u u 


V 

k 

h 3 h 2 

h 1 

h° 

Inv. 

d] 

free 

N f 

free 

d 2 

next 

N 

next 

Y (dB) 

1 

1 

- 

1 

3 

~ — " — ■ — xr 

22.5 

0.890 

8 

- 

- 

1.82 

2 

1 

- 

2 

5 

22.5° 

1.172 

4 

1.476 

32 

3.01 

3 

2 

- 04 

02 

11 

0 

90 

1.476 

16 

- 

- 

4.01 

4 

2 

14 

06 

23 

O 

45 

1.757 

8 

- 

- 

4.77 

5 

2 

- 30 

16 

41 

O 

45 

1.781 

16 

- 

- 

4.83 

6 

2 

- 044 

016 

107 

O 

45 

2.0 

4 

2.085 

48 

5.33 

7 

3 

110 044 016 

317 

O 

45 

2.085 

25 

- 

- 

5.51 


y = 0 dB 
'8 
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Table 21(a): Trellis Coded 3xl6PSK. 


R = 3.67 bit/T, q=0, d 2 = 0.304, N = 12 (3xl6PSK I). 

eff M u u 


V 

k 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N f 

free 

d 2 

next 

N 

next 

Y (dB) 

sig. 

set 

1 

1 

- 

- 

1 

3 

22.5° 

0.304 

4 

- 

- 

0.00 

II 

2 

1 

- 

- 

2 

5 

22.5° 

0.457 

16 

- 

- 

1.76 

II 

3 

2 

- 

04 

02 

11 

22. 5°^ 

0.586 

6 

0.609 

16 

2.84 

I 

4 

3 

14 

04 

02 

21 

45 

0.609 

12 

- 

- 

3.01 

I 


3 

10 

04 

02 

21 

0 

90 

0.609 

8 

- 

- 

3.01 

I 

5 

3 

30 

14 

02 

53 

0 

45 

0.761 

48 

- 

- 

3.98 

I 

6 

3 

050 022 006 

103 

O 

45 

0.890 

12 


- 

4.66 

I 

7 

3 

056 

112 004 225 

0 

45 

0.913 

84 

- 

- 

4.77 

I 


3 

100 050 022 255 

0 

90 

0.913 

76 

- 

- 

4.77 

I 


y = 0 dB 
'8 


Table 21(b): Trellis Coded 3xl6PSK. 


R = 3.33 bit/T, q=l, d 2 = 0.457, N = 8 (3xl6PSK n). 

eff u u 


V 

k 

h 3 

h 2 

h 1 

h° 

Inv. 

d f 

free 

N f 

free 

d 2 

next 

N 

n e xt 

Y (dB) 

sig. 

set 

1 

1 

- 

- 

1 

3 

U 

45 

0.586 

6 

0.609 

16 

1.08 

II 

2 

2 

- 

3 

1 

7 

0 

45 

0.738 

6 

- 

- 

2.08 

II 

3 

2 

- 

06 

02 

11 

O 

45 

0.913 

16 

- 

- 

3.01 

II 


2 

- 

04 

02 

11 

O 

90 

1.043 

24 

- 

- 

3.58 

II 

4 

3 

10 

04 

06 

21 

22. 5° 

1.043 

12 

- 

- 

3.58 

III 


2 

- 

14 

02 

27 

O 

45 

1.172 

12 

1.195 

24 

4.09 

II 

5 

3 

34 

16 

06 

41 

22.5° 

1.172 

4 

- 

- 

4.09 

III 

6 

3 

032 046 006 

103 

22.5° 

1.218 

8 

- 

- 

4.26 

III 

7 

3 

014 

102 044 203 

22.5° 

1.370 

32 

- 

- 

4.77 

III 


3 

006 072 062 

223 

O 

45 

1.476 

8 

- 

- 

5.09 

III 


Y = -1.97 dB 
'8 
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Table 21(c): Trellis Coded 3xl6PSK. 


R = 3.00 bit/T, q=2, d 2 = 0.586, N = 2 (lx8PSK). 

eff u u 


V 

k 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N f 

free 

d 2 

next 

N 

next 

Y (dB) 

sig. 

set 

1 

1 

- 

- 

1 

3 

D 

90 

1.043 

24 

- 

- 

2.50 

II 

2 

1 

- 

- 

2 

5 

0 

90 

1.172 

12 

1.628 

144 

3.01 

II 

3 

2 

- 

04 

02 

11 

22.5° 

1.172 

4 

- 

- 

3.01 

III 

4 

2 

- 

12 

04 

27 

22.5° 

1.628 

32 

- 

- 

4.44 

III 

5 

2 

- 

14 

02 

41 

22.5° 

1.628 

16 

- 

- 

4.44 

III 


2 

- 

22 

14 

43 

0 

45 

1.757 

16 

- 

- 

4.77 

III 

6 

2 

- 

054 

020 

115 

22.5° 

1.757 

8 

2.085 

48 

4.77 

III 


3 

020 

004 

012 

101 

0 

45 

i 

2.0 

6 

2.085 

72 

5.33 

II 


3 

050 

030 

026 

101 

0 ' 

90 

j 

2.0 

6 

2.085 

60 

5.33 

II 

7 

3 

060 

106 

050 

213 

! O ! 

45 

2.0 

6 

2.214 

56 

5.33 

III 


3 

016 

110 

052 

203 

0 

90 

2.0 

6 

2.343 

64 

5.33 

III 


y = 0 dB 

'g 
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Table 22(a): Trellis Coded 4xl6PSK. 


R = 3.75 bit/T, q=0, d 2 = 0.304, N = 24 (4xl6PSK). 

cff u u 


V 

k 

h 4 h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N 

free 

d 2 

next 

N 

next 

Y (dB) 

1 

1 

- 

- 

1 

3 

XT 

22.5 

0.304 

8 

0.457 

64 

0.00 

2 

2 

- 

2 

1 

5 

22.5° 

0.457 

48 

- 

- 

1.76 

3 

2 

- 

04 

02 

11 

22.5° 

0.586 

8 

0.609 

64 

2.84 

4 

3 

- 10 

04 

02 

21 

22.5° 

0.609 

40 

- 

- 

3.01 

5 

3 

- 30 

14 

02 

41 

22.5° 

0.609 

8 

0.761 

288 

3.01 

6 

4 

030 020 052 014 

101 

22.5° 

0.761 

136 

- 

- 

3.98 


y = -1.87 dB 
'8 


Table 22(b): Trellis Coded 4xl6PSK. 


R = 3.50 bit/T, q=l, d 2 = 0.304, N = 4 (2xl6PSK). 

eff 1 u u 


V 

£ 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N r 

free 

d 2 

next 

N 

next 

Y (dB) 

1 

1 

- 

- 

1 

3 

U 

22.5 

0.586 

8 

0.609 

64 

2.84 

2 

2 

- 

2 

1 

5 

22.5° 

0.609 

40 

- 

- 

3.01 

3 

2 

- 

04 

02 

11 

22.5° 

0.609 

8 

0.890 

32 

3.01 

4 

3 

14 

04 

02 

21 

22.5° 

0.890 

16 

- 

- 

4.66 

5 

3 

24 

14 

02 

41 

22.5° 

0.913 

64 

- 

- 

4.77 

6 

3 

014 024 042 

103 

22.5° 

1.172 

24 

1.218 

1088 

5.85 


y = -2.17 dB 
'8 
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Table 22(c): Trellis Coded 4xl6PSK. 


R = 3.25 bit/T, q=2, d 2 = 0.586, N = 8 (4xl6PSK). 

eff u u 



V 

k 

h 4 

h 3 

h 2 

h 1 

h° 

Inv. 

d 2 

free 

N 

free 

d 2 

next 

N 

next 

Y (dB) 


1 

1 

- 

- 

- 

1 

3 

U 

22.5 

0.609 

8 

0.890 

32 

0.17 


2 

2 

- 

- 

3 

1 

5 

22.5° 

0.890 

16 

- 

- 

1.82 


3 

2 

- 

- 

06 

02 

11 

22.5° 

1.172 

24 

1.195 

64 

3.01 



2 

- 

- 

02 

06 

11 

22.5° 

1.172 

24 

1.218 

64 

3.01 


4 

3 

- 

04 

06 

12 

21 

22.5° 

1.172 

8 

1.218 

32 

3.01 


5 

4 

10 

04 

06 

22 

41 

22.5° 

1.218 

16 

- 

- 

3.18 

— 

6 

4 

050 030 024 016 

101 

22.5° 

1.499 

72 

- 

- 

4.08 


y g = 0.35 dB 


Table 22(d): Trellis Coded 4xl6PSK. 


R = 3.00 bit/T, q=3, d 2 = 0.586, N = 2 (lx8PSK). 

eff u u 


V 

E 

h 3 

h 2 

h 1 

h° 

Inv. 

d f 

free 

N f 

free 

d 2 

next 

N 

next 

Y (dB) 

1 

l 

- 

- 

1 

3 

o — 

45 

1.172 

24 

1.218 

64 

3.01 

2 

2 

- 

2 

3 

5 

22.5° 

1.172 

8 

1.218 

32 

3.01 

3 

3 

02 

04 

03 

11 

22.5° 

1.218 

16 

- 

- 

3.18 

4 

3 

04 

10 

06 

21 

22.5° 

1.781 

48 

- 

- 

4.83 

5 

3 

22 

16 

06 

41 

22.5° 

1.804 

24 

- 

- 

4.88 


3 

24 

14 

02 

43 

O 

45 

1.827 

64 

- 

- 

4.94 

6 

3 

050 024 006 

103 

22.5° 

2.0 

8 

2.343 

64 

5.33 


y 8 = o dB 
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1 Introduction 


Demands for higher information rates (spectral efficiencies) in satellite com- 
munications are such that this traditionally power limited channelis now 
considered to be bandlimited as well. This has made high speed satellite 
channels a natural candidate for the application of Trellis Coded Modulation 
(TCM). Introduced by Ungerboeck [1] ,[2], TCM achieves significant asymp- 
totic coding gains relative to an uncoded system of the same information 
rate with no bandwidth expansion. Due to the extensive use of nonlinear 
amplifiers, such as traveling wave tubes, in satellite communications systems 
only TCM schemes using MPSK signal constellations are appropriate. A 
class of multidimensional (multi-D) MPSK trellis codes has been proposed 
for use in NASA high speed satellite communication links [3]. These codes 
have greater flexibility in achievable information rate, better rotational in- 
variance properties, and higher speed decoder implementations than regular 
two dimensional (2-D) Ungerboeck codes. 

A typical satellite communication system consists of a Reed-Solomon(RS) 
outer code, a symbol interleaver, and an inner code (Figure 1). The satellite 
channel is usually well modeled by the zero mean Additive White Gaussian 
Noise (AWGN) process with two sided power spectral density N a /2. In one 
mode of operation, the RS code and interleaver are bypassed and the in- 
formation bits are directly encoded by the inner code and transmitted. In 
this case, a bit error rate (BER) of ~ 10 -5 out of the inner decoder is con- 
sidered reasonable. Under poor channel conditions or when extremely low 
error rates are required, the entire concatenated coding system is used. For 
a given ( N , K) RS code over GF{ 2 b ) and assuming ideal interleaving, the 
performance of the system is determined by the symbol error rate (SER) out 
of the inner decoder, where the symbol size, b, is determined by the RS code. 
Thus, the choice of the inner code is critical when designing such a system. 

In this paper, the performance of multi-D MPSK trellis codes as the 
inner code in a satellite communication system is studied through extensive 
computer simulation. In Section 2, a brief review of the results in [3] is 
presented in order to introduce multi-D MPSK signal sets and the multi-D 
encoder structure. In Section 3, the BER performance of multi-D MPSK 
trellis codes on the AWGN channel with soft decision Viterbi decoding is 
studied. This section includes an example of the calculation of the Zehavi and 
Wolf transfer function bound and results on the minimum truncation length 
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required for Viterbi decoding of MPSK trellis codes. Finally, the performance 
of MPSK trellis codes in a concatenated coding system is studied in Section 
4 through the use of computer simulation to obtain the SER out of the inner 
decoder. 


2 Multi-Dimensional Trellis Codes 

The fundamental concept of 2-D TCM is that by trellis coding onto an ex- 
panded signal set relative to that needed for uncoded transmission, a larger 
free distance and an associated improvement in performance may be achieved 
without bandwidth expansion or a reduction in information rate. In [1] 
and [2], Ungerboeck considered the 2-D MPSK case where uncoded 2 m PSK 
was expanded into a rate m/(m + 1), 2 m+1 PSK trellis code. In multi-D 
MPSK TCM, the same principal is applied to a signal constellation, denoted 
LxMPSK, constructed from the Cartesian product of L 2D-MPSK signal sets. 
Each branch of the trellis is then labelled with an ordered sequence of L 2-D 
signals and a 2L-dimensional signal is realized by transmitting a sequence 
of L 2-D signals. The efficient labelling of trellis branches is accomplished 
through a systematic binary set partitioning technique based on block codes 
and their cosets. 

The class of LxMPSK codes studied here has the same schematic encoder 
diagram introduced by Forney et.al.[4] for Ungerboeck codes with one slight 
modification. As shown in Figure 2, the multi-D encoder allows the q least 
significant bits to be fixed, nominally to 0, during each coding interval. Thus, 
whereas 2D-MPSK TCM allowed only rate m/m + 1 codes with information 
rates ( R ) of 1.0, 2.0, and 3.0 bits per modulation interval(T) for QPSK, 
8PSK, and 16PSK, respectively, multi-D MPSK allows the construction of 
rate (Xxlog 2 (M — q — l))/(Ioclog 2 M) codes with a variety of information 
rates ranging from 1.0 bit/T to 3.75 bits/T. 

As an example, consider the familiar 4-state Ungerboeck 8PSK code 
shown in Figure 3 with information rate 2.0 bits/T. The fundamental trellis 
structure of this code is determined by the rate 1/2 convolutional encoder 
and is shown in Figure 4. During each coding interval, the trellis encoder 
receives m = 2 information bits, of which m = 1 are encoded by the rate 
1/2 convolutional encoder, and maps them to a single 8PSK signal. This 
results in the trellis structure of Figure 4 with each transition replaced by 
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2 m -m _ 2 1 parallel transitions. The same rate 1/2 convolutional encoder 
can be used to construct a 2x8PSK code with information rate 2.5 bits/T. 
In this case, the trellis enco der receives m — 5 information bits per coding 
interval, of which the least significant bit is encoded by the convolutional 
encoder. The m + 1 = 6 bits are then mapped to a single 2x8PSK signal. 
This results in a rate 5/6 code with the trellis structure of Figure 4 where 
each branch is replaced with 2 m_ * = 2 4 = 16 parallel transitions and each 
transition is labelled with a 2x8PSK signal. If the first q = 1 bits into the 
trellis encoder are now set to 0, a rate 4/6, 2x8PSK code with information 
rate 2 bits/T is obtained. This code has the trellis structure of Figure 4 with 
2 4-1 = 8 parallel transitions and each transition is labelled with a 2x8PSK 
signal from the left branch of the first level of the 2x8PSK binary partition 
(see Pietrobon [3] et. al. for a discussion of multi-D signal set partitioning). 
It should be noted that, though this example illustrates the concept of these 
LxMPSK trellis codes, the best LxMPSK codes will in general not have the 
same convolutional encoder as a 2-D Ungerboeck code. 


3 Code Performance and Simulation Results 


Having introduced the basic concepts of LxMPSK trellis codes, it is now of 
interest to investigate their performance. It is common practice to measure 
the performance of trellis codes in terms of the asymptotic coding gain rela- 
tive to an uncoded system with the same rate. In this case, the coding gain 
is given by 


7=10 log 


10 


Al 

~freercM 


^/rcCuncoded J 
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The asymptotic coding gain, however, is accurate only as the SNR becomes 
very large and does not provide much insight into a codes real coding gain 
at a particular SNR or error rate. In order to determine real coding gains, 
analytical bounds or simulations must be used. 

For linear convolutional codes, the most common analytical method for 
determining code performance is the transfer function bound. This method 
uses a modified state diagram to enumerate the complete set of possible error 
paths from the correct sequence. Inherent in this technique is the assumption 
that considering the all zero sequence as the correct sequence entails no loss 
of generality. For most TCM schemes, the mapping from code sequences to 
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channel symbols is non-linear and the set of error paths is dependent on which 
code sequence was sent. Thus, it is necessary to enumerate the set of all error 
paths for each possible correct code sequence. This requires a state diagram 
using 2 2l/ states and quickly becomes Impractical[5]. Zehavi and Wolf [6] 
have shown that for a large class of trellis codes, including those discussed 
in [1],[2], and [3], a transfer function upper bound on the average probability 
of error can be calculated using a modified state diagram having only 2" 
states. As with convolutional codes, this procedure becomes intractable as 
the number of state increases. For multi- D trellis codes, the Zehavi and Wolf 
method becomes even more complex due to an exponential increase in the 
number of error vectors whose distance profile from each possible code vector 
must be enumerated. The calculation of the Zehavi and Wolf bound and its 
computational difficulties is best manifest through an example. 

Consider, the 8PSK, // = 2, Ungerboeck code introduced in Section 2 
and shown in Figure 3. A necessary condition for the application of this 
bound is that both subsets at the first level of the binary partition tree must 
have the same distance profile with respect to a binary error vector E. In 
this case, it is simple to check that the two subsets A={S 0 , S 2 , S 4 , S 6 } and 
k={S u S 3 ,S 5 ,S 7 } have the same profile for all 8 possible three bit error 
vectors. For example, the error vector Oil has squared Eulcidean distances 
{3.4142,0.5858,3.4142,0.5858} from A and {0.5858,3.4142,0.5858,3.4142} 
from A. With this condition, the distance profile for each of the 8 error 
vectors with respect to either subset can be determined as shown in Table 
1. Note, that for a LxMPSK code this involves finding the distance profiles 
of error vectors with respect to M L / 2 signal points. From Table 1, 

the distance polynomial of each error vector can be obtained as in Table 2, 
where the power of W is the squared Euclidean distance and the coefficient 
is the number of error paths with that distance. This information is used 
to construct a modified state diagram that enumerates the multiplicity and 
distance of all possible error paths. 

The state diagram of the R=l/2 convolutional code used in this trellis 
code is shown in Figure 5. The diagram of Figure 5 is transformed into 
modified state diagram that enumerates all possible error paths by labelling 
each branch with a polynomial consisting of terms of the form L x PW a 
where i is the number of symbols, j is the number of nonzero informa- 
tion bits, and a is the squared Euclidean distance. The key to labelling 
the branches of the modified state diagram is to remember that each tran- 
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sition in the convolutional encoder’s state diagram represents a group of 
parallel transtions for the trellis code. Tor example, when considering the 
transiton from state 1 to state 3 in Figure 5 the polynomials of the error 
vectors Oil and 111 must be taken into account yielding the polynomial 
2LI(2W° 58579 + 2 IF 3 - 4142 ) + 2LI 2 (2W 0 58579 + 2W 3-4142 ) as the appropriate 
label. Repeating this procedure for each branch gives the modified state 
diagram of Figure 6. . „ . 

Regarding this as a signal flow graph, the transfer function is found to be 


T{W,LJ) = iLIW A + 


L 4 I 4 a 2 p 2 + L 3 I 2 a 2 j(l - I 7 ) 

1 - L 3 I 2 0 2 8 - L 2 1 8 - I 7 + L 3 P8 


(2) 


where a = 4(7 + 1 )W 2 , j3 = 2(7 + 1 )(IF °- 58579 + W 3A142 ), 7 = 4 (IF 058579 + 
/VF 3-4142 ), and 8 = 4(1 + IW A ). Using this transfer function, Zehavi and 
Wolf show that the bit error rate for a trellis code with maximum likelihood 
decoding (MLD)is upper bounded by 
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where d 2 is the normalized squared free Euclidean distance of the code. The 
Zehavi and Wolf bound and corresponding simulation results are plotted in 
Figures 7 and 8 for the Ungerboeck code discussed above and a 2x8PSK, 
R=2.5 bits/T, 4-D code from [3], respectively. As expected, the bound and 
the simulation results conform closely except at low SNR’s where the union 
bound becomes relatively loose. This of practical interest since many satellite 
communications systems operate at low SNR’s and concatenated systems are 
sensitive to small deviations in the performance of the inner code. 

Simulation results for these codes were obtained via a stochastic or Monte 
Carlo simulation of a soft decision MLD Viterbi decoder. This involves sim- 
ulating the AWGN channel with a computer generated pseudo-random num- 
ber generator that approximates the Gaussian distribution. The Marsaglia- 
Bray[7] method was used to implement the Gaussian random number gener- 


ator from the available uniform random number generators. Because trellis 
codes are non-linear, a number of nonzero information sequences are con* 
sidered in each simulation and the code performance is determined as the 
average over these sequences. 
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Figure 9 shows the simulation results for a group of v = 2, Lx8PSK 
trellis codes with information rate 2.0 bits/T and <Pj ree = 4.0. These codes 
show a coding gain from 2.0dB to 2.4dB at ~ 10 -5 compared to uncoded 
QPSK. The 4, 6, and 8 dimensional codes suffer a 0.2dB to 0.4dB loss in 
performance compared to the Ungerboeck 8PSK code due to higher path 
multiplicities and denser distance spectra. For example, the 4x8PSK code 
has a multiplicity N j Tee = 12 of the minimum free distance path (dy ree = 4.0) 
and a second spectral line at d\ = 4.686 with multiplicity iV 2 = 32. On the 
other hand, it is possible to implement a 4x8PSK decoder at 4 times the speed 
of the Ungerboeck code by using read only memories (ROM’s) to decode 
the parallel transitions. Analagous simulation results were obtained for the 
multi-D Lx4PSK codes with information rate 1 bit/T and the Lxl6PSK 
codes with information rate 3 bits/T. Figure 10 shows simulation results for 
v = 4 and v = 6, R = 2.67 bits/T, 3x8PSK codes compared to a v — 4 8PSK 
periodically time varying trellis code (PTVTC) [8] of the same effective rate. 
The 3x8PSK trellis codes show an improvement of 0.5 dB for v — 4 and 
1.0 dB for v = 6 and are capable of a decoder implementation that would 
operate at three times the speed of the PTVTC decoder. Simulation results 
for a number of the codes given in [3] are contained in Appendix I . 

The ideal Viterbi decoder retains the surviving path into each encoder 
state and waits until the end of the current code sequence before decoding 
any information bits. In actual implementations, finite memory limitations 
allow the decoder to retain only the past r branches of information bits for 
each state. The effect of this constraint, called the truncation length of the 
decoder, is that additional decoding errors are introduced other than those 
that would occur in an ideal MLD Viterbi decoder. For a particular code, 
the minimum truncation length is the minimum value of r for which the per- 
formance of the truncated Viterbi decoder will asymptotically approach that 
of the ideal MLD decoder. In [9] the Zehavi and Wolf bound is modified to 
account for the truncation length of the Viterbi decoder. Using the modified 
bound, the performance of a code for a particular r can be calculated as well 
as the minimum truncation length required for asymptotically ML perfor- 
mance. In Figure 11, the modified bound and simulation results are plotted 
for varying r for the R=2 bits/T, v = 2, 8PSK Ungerboeck code. Figure 
12 shows the truncation length simulation results for the 3x8PSK, R=2.33 
bits/T, v = 4 code currently being implemented. In [9], the minimum trun- 
cation length of this code was calculated to be 11 branches. 
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4 Concatenated Coding and LxMPSK Codes 

Many satellite communications links use concatenated coding in order to 
achieve extremely low BER’s over noisy channels with reasonable code com- 
plexity. The class of Reed-Solomon (RS) block codes are often chosen as the 
outer code in these systems due to their optimal distance properties, high 
rates, and the availability of relatively simple decoding algorithms. Addi- 
tionally, an interleaver is usually placed between the inner and outer code 
to ensure that errors in the RS symbols are independent. For a concate- ■ 
nated coding system using a t-error-correcting ( N,K ) RS outer code and 
ideal interleaving, the overall bit error rate is well approximated by 


Pb 


2t -f" 1 
N 




(4) 


where p, is the RS symbol error rate out of the inner decoder. Computer 
simulation is necessary to get an accurate measure of p» for a given inner 
code. 

When using convolutional inner codes, it was shown in [10] that choosing 
an RS code over GF( 2 6 ) matched to a byte oriented inner convolutional codes 
of rate k/b resulted in an improvement in the overall sytem BER by nearly 
an order of magnitude compared to unmatched inner codes of the same rate. 
Costello and Deng[ll] have suggested that when using multi-D LxMPSK 
inner codes and RS outer codes similar results might apply. In Figure 13, 
the BER and symbol error rate (SER) for a v = 2, R = 2.0 bits/T, 4x8PSK 
code and the comparable 2-D Ungerboeck code are shown for a symbol size 
of 8 bits, i.e. the RS code is over GF( 2 8 ). The results indicate that the BER 
penalty of 0.2dB to 0.4dB that multi-D codes suffer compared to 2D codes, as 
noted above, disappears when SER is considered. Figure 14 shows the SER 
simulation results for the 3x8PSK, v = 4, R - 2.33 bits/T code currently 
being implemented versus an uncoded 3x8PSK system with the same rate. 
Thus, the multi-D codes appear to be better suited for concatenated coding 
than the 2-D codes. 
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5 Conclusion 


The performance of a class of multi-D LxMPSK codes was studied through 
extensive computer simulation. When compared to 2D Ungerboeck codes 
with the same information rate, the multi-dimensional codes performed slightly 
worse due to their high path multiplicities and dense spectra. Additionally, 
they have a higher trellis complexity due to the larger number of parallel 
transitions. However, these codes have better rotational invariance proper- 
ties and are capable of achieving much higher decoder speeds. Furthermore, 
multi-D LxMPSK codes can be constructed using rates not possible with 
conventional 2D-MPSK codes. It was also shown that the byte oriented 
structure of the multi-D codes compensated for the slight loss in BER so 
that a concatenated system with the proper RS outer code would perform 
equally as well as a 2D code. 
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Table 1: Distance profile for the 8-PSK, v = 2 Ungerboeck code. 
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Table 2: Branch Weights for the 8-PSK, u = 2, Ungerboeck code. 


E = ( E 2 E l E ° ) 

Weight Profile 

000 

4 

001 

4^0.58579 

010 

4W 2 

Oil 

2py0. 58579 2^3-4142 

100 

4 IV 4 

101 

4^3.4142 

110 

4W 2 

111 

2^0.58579 _j_ 2py3.4142 


12 



Figure 1: Block Diagram of a Typical Satellite Communications System 
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Figure 2: Schematic Representation of a Multi-D Trellis Encoder 













Figure 3: Encoder for the 4 State, 8PSK, R=2/3 Ungerboeck Code 
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Figure 4: Trellis Diagram for the Rate 1/2 Convolutional 
Encoder of Figure 3. 
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Figure 8: The Performance of a 4 State, 2x8PSK, R=2.5 bits/T Multi-D Trellis Code 
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Figure 10: Simulation Results for Trellis Codes with R=2.67 bits/T. 
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Figure 11: Truncation Length Results for the 4 State, 8PSK, Ungerboeck Code. 
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Figure 14: Symbol Error Rate Simulation Results for a RS(255,223) Outer Code. 






7 Appendix I 

This appendix contains the simulation results for a number of LxMPSK trellis 
codes. All table numbers refer to reference [3]. 
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Figure A.4: 3x4PSK, R=1.67 bits/T Simulation Results (Table 13a) 
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Figure A.9: 4x4PSK, R=1.25 bits/T Si 
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Figure A.10: 4x4PSK, R=1.00 bits/T Sin 
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Figure A. 11: 8PSK, R=2.00 bits/T Simulation Results (Table 15) 
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Figure A.12: 2x8PSK, R=2.50 bits/T Simulation Results (Table 16a) 
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Figure A. 13: 2x8PSK, R=2.00 bits/T Simulation Results (Table 16b) 
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Figure A.14: 3x8PSK, R=2.67 bits/T Simulation Results (Table 17a) 
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Figure A. 15: 3x8PSK, R=2.33 bits/T Simulation Results (Table 17b) 
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Figure A. 17: 4x8PSK, R=2.75 bits/T Simulation Results (Table 18a) 
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Figure A. 18: 4x8PSK, R=2.50 bits/T Simulation Results (Table 18b) 
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Figure A. 19: 4x8PSK, R=2.25 bits/T Simulation Results (Table 18c) 
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Figure A.24: 3xl6PSK, R=3.6 








I- 01 



3*^ jojih 


Figure A.25: 3xl6PSK, R=3.33 bits/T Simulation Results (Table 21b) 









